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SUMMARY 


-  An  improved  method  is  presented  for  calculating  rotor  system  over¬ 
all  vortex  noise  and  frequency  spectra  for  stalled  and  unstalled  rotors. 
Correlation  of  measured  and  predicted  vortex  noise  was  evaluated  using 
two  rotor  systems  operating  over  a  wide  range  of  speeds  and  thrusts. 
Correlation  was  found  to  be  excellent.  Blade  tip  planform  studies 
revealed  significant  vortex  noise  reductions  with  tapered  tips. 

A  new  procedure  is  also  derived  for  calculating  near  and  far  field 
rotor  rotational  noise  with  nonuniform  inflow.  The  method  extends  the 
standard  steady  load  method  by  including  the  effects  of  harmonic  airloads. 
Correlation  studies  were  conducted  using  an  H-34  helicopter.  Agreement 
between  low  frequency  measured  and  predicted  noise  was  good.  However, 
correlation  with  high  harmonic  rotational  noise  was  poor.  This  is  prob¬ 
ably  due  to  inadequate  definition  of  high  harmonic  airloads. 

Presented  results  establish  the  importance  of  high  harmonic  rota¬ 
tional  noise  for  detectability  and  loudness,  and  further  work  is  recom¬ 
mended  to  more  accurately  define  high  harmonic  blade  loading.  Since 
an  airload  measurement  program  is  being  conducted  on  the  NH-3A,it  is 
recommended  that  a  correlation  program  be  conducted  to  more  fully 
evaluate  the  accuracy  of  die  presented  noise  analysis  program  using  the 
NH-3A  airload  results. 

This  study  was  performed  for  single  rotor  systems  only,  and  in  its 
present  form  is  not  directly  applicable  to  systems  with  multiple  rotors  in 
juxtaposition. 
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PHASE  I 


MAIN  ROTOR  VORTEX  NOISE  DURING  UNIFORM  INFLOW 
INTRODUCTION 

The  task  of  Phase  I  was  to  develop  an  improved  procedure  for  pre¬ 
dicting  single  rotor  helicopter  main  rotor  vortex  noise  under  conditions 
of  uniform  inflow.  This  task  was  accomplished  by  a  review  of  vortex 
noise  theory,  development  of  a  formula  and  a  program  of  rotor  stand 
measurements.  Although  the  proposed  work  covered  only  the  testing  of 
one  main  rotor  system  under  three  conditions  each  of  disk  loading  and 
rotor  tip  speed,  the  availability  of  another  rotor  system  provided  an 
opportunity  to  check  analytical  accuracy. 

Previous  studies  of  vortex  noise  generation  by  helicopter  rotor  sys¬ 
tems  have  been  limited  to  prediction  of  the  overall  levels  for  systems 
operating  out  of  stall.  Very  little  has  been  reported  on  prediction  of  the 
frequency  spectrum  for  conditions  either  in  or  below  stall.  However, 
the  importance  of  the  spectral  distribution  of  the  noise  in  assessing  the 
effects  of  rotor  design  on  such  psycho-acoustic  factors  as  detectability 
and  hearing  damage  requires  the  estimation  of  the  spectrum  shape. 
Consequently,  the  individual  spectra  were  studied  in  detail. 

This  report  presents  the  results  of  this  study  including  an  improved 
method  of  calculation  of  the  noise  levels  and  normalized  spectra  for 
operation  in  and  out  of  stall.  Sample  spectra  are  given  for  both  condi¬ 
tions.  The  improved  method  consists  of  corrections  to  Harvey  Hubbard's 
formula  (Reference  16)  to  account  for  changes  in  lift  coefficient.  The 
lift  coefficient  is  now  a  variable,  and  the  actual  blade  area  is  used  rather 
than  a  calculated 'effective  blade  area.  The  coefficients  to  be  used  in 
this  formula  have  been  established  empirically  from  CH-3C  and  CH-53A 
rotor  test  stand  data. 

VORTEX  NOISE  GENERATION 


The  subject  of  sound  generation  by  fluid  flow  has  been  studied  inten¬ 
sively  by  a  number  of  investigators.  Recent  advances  in  the  understand¬ 
ing  of  how  sound  is  generated  and  propagated  by  unsteady  aerodynamic 
phenomena  such  as  vortices  have  permitted  rigorous  mathematical  treat¬ 
ment.  Reference  26  contains  a  complete  physical  and  mathematical  des¬ 
cription  of  the  phenomena  of  vortex  sound. 
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Briefly,  vortex  sound  is  generated  by  the  fluid  force  on  an  object 
arising  from  the  formation  and  shedding  of  vortices  in  the  flow  past  it. 
This  results  in  a  dipole  form  of  radiation  in  which  the  strength  of  the 
source  is  proportional  to  the  sixth  power  of  the  free- stream  velocity. 

For  a  rod  or  bluff  body  operating  at  low  Reynolds  numbers,  the  vor¬ 
tices  are  shed  alternately  from  each  side  of  the  rod  in  regular  vortex 
shedding.  The  shedding  of  vortices  causes  fluctuating  lift  and  drag.  The 
sound  associated  with  the  fluctuating  drag  is  much  weaker  and  is  double 
the  frequency  of  that  associated  with  lift. 

In  the  usual  range  of  Reynolds  numbers  (approximately  5  x  10 6) 
associated  with  a  rotating  wing,  the  sound  frequency  is  given  by 

«) 

where  V  is  the  free-stream  velocity  and  d  is  the  thickness  of  the  rod  or 
bluff  body.  The  Strouhal  number,  $4,  has  been  experimentally  deter¬ 
mined  to  have  a  value  of  0. 28  (Reference  28).  Since  there  is  a  different 
velocity  associated  with  each  station  over  the  span,  there  is  a  broad 
band  of  frequencies.  However,  the  intensity  of  sound  is  proportional  to 
the  sixth  power  of  the  velocity  V,  so  that  the  frequencies  of  interest  are 
associated  with  the  area  near  the  tip  where  the  velocity  is  highest. 

Table  I  shows  the  calculated  frequencies  and  their  relative  intensity  for 
the  CH-53A  rotor  blade.  From  Table  I,  it  can  be  seen  that  the  most 
intense  sound  appears  in  the  octave  from  300  to  600  cps.  For  the  pur¬ 
poses  of  this  report,  all  of  the  sound  from  150  to  9600  cps  was  con¬ 
sidered  to  be  vortex  sound. 

Although  vortex  sound  is  the  principal  source  of  medium  and  high 
frequency  noise,  boundary  layer  turbulence  and  noise  from  blade  irregu¬ 
larities  can  contribute  to  the  overall  noise.  Turbulence  on  the  blade  is 
another  source  of  frequency  broadening.  The  turbulence  causes  an 
irregular  rate  of  shedding.  If  the  blade  is  operating  in  a  region  beyond 
the  onset  of  tip  stall,  there  is  a  sharp  rise  in  the  1200-  to  2400-cps  octave 
band  associated  with  the  flow  separation  in  the  region  near  the  tip.  This 
flow  separation  can  be  associated  with  blade  slap  as  shown  in  phase  III. 

CALCULATION  OF  VORTEX  NOISE 


Overall  Level 

A  theory  of  vortex  sound  was  first  presented  by  E.  Y.  Yudin  (Refer¬ 
ence  33).  Yudin's  theory  for  the  radiation  of  vortex  sound  from  rotating 
rods  was  based  on  a  dimensional  analysis  of  the  flow  parameters  around 
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the  rods.  In  developing  his  theory,  Yudin  assumed  that  a  rigid  body  in  a 
moving  fluid  has  forces  impressed  upon  it  due  to  the  shedding  of  vor¬ 
tices.  The  sound  radiation  is  the  same  in  this  case  as  that  from  equal 
but  opposite  forces  acting  directly  on  an  otherwise  still  fluid. 

Curie  (Reference  4),  in  his  analysis  of  the  influence  of  solid  bounda¬ 
ries  on  aerodynamic  sound,  showed  that  Yudin's  result  could  be  obtained 
formally  from  Lighthill's  theory  (Reference  19)  in  which  the  total  fluid 
stress  is  associated  with  the  local  dipole  strength.  Curie  has  thus 
shown  that  Yudin's  relationship  was  correct.  Yudin's  result  can  be 
written  as  follows: 

W# s const. ^3T  (C%  Sjf  Vp  r\\  (2) 

In  using  equation  2  to  calculate  vortex  noise  from  rotating  wings, 
some  of  the  terms  can  be  considered  to  be  constant  because  their  range 
of  permissible  values  is  too  small  to  affect  the  results  appreciably. 

These  are  p,  C  ,  C€  and  S*.  Although  ,  die  form  drag  coefficient, 
appears  in  the  equation  to 'the  fourth  power,  the  range  in  lift  between  an 
unloaded  aircraft  and  a  fully  loaded  aircraft  is  small  enough  to  accept 
an  average  value  for  Q,  for  approximate  calculations.  Lumping  the 
constant  terms  then  gives  the  following: 

W*  5  const  Vp  r  If  (3) 

The  two  helicopter  rotor  systems  tested  had  untapered  0012  airfoils. 
Therefore,  a  constant  proportion  exists  between  the  chord  and  the  thick¬ 
ness  n  .  Thus,  the  sound  power  is  proportional  to  the  chord  length  times 
the  blade  length  r,  or  blade  area.  Furthermore,  since  the  noise  from 
each  blade  is  additive,  the  total  sound  power  will  be  proportional  to  die 
blade  area  of  the  rotor  system.  This  simplifies  the  equation  to  the 
fol’  owing: 

)Nk  s  const.  Ag  Vp  (4) 

Equation  4  can  new  be  compared  directly  to  Hubbard's  formula 
(Reference  16)  which  was  used  to  calculate  the  sound  pressure  level  of 
the  vortex  noise  from  propellers,  which  is  as  follows: 


Hubbard's  measurements  were  made  at  a  distance  of  300  feet  and  an 
angle  of  105*  from  die  axis  of  rotation.  The  constant  K  in  equation  5  was 
determined  empirically  for  the  measurement  position.  The  terms 
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10  log  and  10"^  give  the  sound  pressure  level  in  decibels. 

Hubbard's  formula  is  based  on  Cj_  =  0.  4.  He  adjusts  the  formula 
for  other  values  of  Cj_  by  using  an  effective  blade  area,  A».  Intensive 
analysis  of  experimental  rotor  test  data  indicated  that  greater  accuracy 
could  be  attained  by  using  the  actual  blade  area  and  coefficient  of  lift. 

The  test  data  yielded  a  value  of  6. 1  x  10'27  for  the  constant  K  of  equation 
5  for  use  in  a  modified  equation.  Variations  in  lift  for  the  modified 
equation  are  accounted  for  by  addition  of  the  term  20  log  C{_/0. 4  where 
the  0. 4  is  the  coefficient  of  lift  used  by  Hubbard  in  equation  5.  The 
resulting  equation  is  as  follows:  g 

SPL*  io  log  *  20 1>9  <*•> 

This  equation  may  be  rewritten  in  a  more  convenient  form  for  sea  level 
70°F  conditions  as  follows: 

5 PL  9  lofelog  Va7  +  Z  IdjT  -  log  Ag  -  j.jrJ  (6b) 

Overall  vortex  noise  levels  were  calculated  for  the  CH-3C  and 
CH-53A  rotor  systems  by  means  of  equations  5  and  6.  The  results  of 
these  calculations  are  shown  in  Table  II.  All  numbers  have  been 
rounded  off  to  the  nearest  db.  As  can  be  seen  in  Table  II,  the  calculated 
levels  agree  with  the  measured  levels  within  2  db.  Also,  it  can  be  seen 
that  equation  5  is  sufficiently  accurate  for  all  but  the  most  stringent 
requirements. 

Since  equations  5  and  6  are  for  a  distance  of  300  feet  from  the  center 
of  rotation,  the  calculated  values  in  Table  II  have  been  increased  by  2. 6 
db  to  correct  for  the  distance  of  225  feet  to  the  point  of  measurement. 
The  usual  distance  corrections  (20  log  300/P)  can  be  made  for  distances 
other  than  300  feet.  * 

Spectrum  Shape 


The  spectrum  shape  of  a  blade  operating  out  of  stall  is  shown  in 
Figure  19.  This  condition  is  present  at  low  angles  of  attack  at  the  tip. 
The  peak  frequency  is  determined  by  the  Strouhal  number  which  is 
defined  as  follows: 

V-w  (7) 

In  the  usual  range  of  Reynolds  numbers  for  a  helicopter  rotor,  the 
Strouhal  number  is  0. 28.  In  keeping  with  Yudin's  "round  rod"  approach, 
*  A  sample  calculation  is  shown  in  Appendix  !. 
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the  projected  blade  thickness  1/  Is  defined  by  the  following  equation: 

h/=  b  cosot  +  a  sin  a  (8) 

where  b  is  the  blade  thickness,  A  the  chord  length  and  a the  angle  of 
attack. 

Tip  Stall 

When  unsteady  aerodynamic  forces  appear  near  the  tip  of  a  blade 
due  to  the  occurrence  of  either  stall  or  drag  divergence,  there  is  a 
definite  change  in  the  shape  of  the  vortex  noise  level  frequency  spectrum. 
A  portion  of  the  rotor  test  data  acquired  was  taken  while  portions  of  the 
blades  were  experiencing  drag  divergence  (and  probably  stall  as  well), 
and  the  general  spectrum  obtained  is  shown  in  Figure  20.  Compare  this 
spectrum  with  that  of  Figure  19  where  no  separated  flow  is  present  and 
it  is  evident  that  the  separated  flow  has  caused  a  rise  in  the  levels  of  the 
octaves  above  the  peak  octave.  Using  these  curves  (Figures  19  and  20), 
one  may  determine  whether  stall  or  drag  divergence  is  present  in  a 
measured  vortex  noise  spectrum.  For  prediction  purposes,  the  proper 
spectrum  may  be  selected  chi  the  basis  of  aerodynamic  stall  and  drag- 
divergence  criteria. 

The  exact  point  where  tip  stall  begins  is  difficult  to  determine  accu¬ 
rately.  At  that  point,  however,  there  is  a  deterioration  of  the  lift /drag  . 
ratio.  Figure  21  is  a  drag-divergence  curve  for  the  CH-3C  blade.  The 
experimental  data  points  plotted  indicated  that  some  of  the  Mach  number- 
angle  of  attack  combinations  are  in  the  region  of  drag  divergence. 

In  Figure  22  the  difference  between  the  300  to  600  octave  levels  and 
the  1200  to  2400  octave  levels  is  plotted  against  the  blade  tip  pitch  for  the 
same  data  points  of  Figure  21.  The  rise  of  the  1200- to  2400-cps  octave 
relative  to  the  300  to  600  octave  corresponds  with  the  data  points  that 
lie  above  the  drag-divergence  curves.  The  rise  appears  less  abrupt  for 
the  high  Mach  numbers. 

EXPERIMENTAL  TEST  PROGRAM 


General 


The  object  of  the  test  program  was  to  measure  the  vortex  noise 
radiated  by  the  rotor  system  over  a  range  of  speeds  and  loads.  Two 
rotor  systems  were  measured  at  three  different  speeds  each  and  at 
several  different  thrusts  at  each  speed.  Test  instrumentation  is 
described  in  Appendix  U. 


/ 
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Description  of  Rotor  Systems 


Both  of  the  two  systems  tested  used  untapered  0012  airfoils  and  had 
square  tips.  The  CH-3C  system  had  five  rotor  blades  and  a  diameter  of 
62  feet.  The  CH-53A  system  had  six  blades  and  a  diameter  of  72  feet. 
Although  the  program  originally  called  for  tests  of  only  one  rotor  system, 
the  opportunity  to  check  the  results  with  another  quite  different  arrange¬ 
ment  was  valuable.  The  physical  characteristics  of  each  system  are 
shown  in  Table  III. 

Acoustical  Measurements 


The  rotors  were  tested  on  a  whirl  stand  approximately  70  feet  above 
the  ground.  Figure  1  shows  the  test  arrangement.  The  ground  plane 
around  the  test  stand  was  covered  with  heavily  matted  dry  grass  except 
for  a  roadway.  The  microphone  was  mounted  on  a  wire  fence  225  feet 
to  the  east  of  the  stand  and  about  six  feet  above  the  ground. 


Figure  1.  Vortex  Noise  Test  Arrangement  and  Rotor  Test  Stand 

There  was  a  grassy  ridge  beyond  the  fence  and  a  wooded  hill  to  the  west 
of  the  stand.  All  measurements  were  made  at  the  same  microphone 
position.  It  is  assumed  that  the  measurements  were  the  direct  field 
from  the  test  rotor  because  (1)  the  portion  of  the  spectrum  of  interest 
was  above  150  cps,  (2)  the  test  rotor  was  70  feet  above  the  ground,  and 
(3)  the  ground  was  covered  with  an  absorbent  grass  cover. 

Since  the  sound  fields  are  strongly  affected  by  winds,  measurements 
were  made  when  the  ambient  wind  was  below  5  knots.  Thus  the  measured 
data  should  closely  approximate  free  field  conditions  unaffected  by  non- 
uniform  inflow. 
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PHASE  II 


HELICOPTER  ROTATIONAL  NOISE  UNDER  CONDITIONS  OF 

NONUNIFORM  INFLOW 


INTRODUCTION 


Rotational  noise,  which  is  comprised  of  discrete  frequencies  at 
multiples  of  che  blade  passage  frequency,  as  opposed  to  vortex  noise, 
which  is  broad-band  random  noise  centered  about  the  Strouhal  frequency, 
was  first  defined  by  Gutin  (Reference  11)  in  1948.  Since  the  establish¬ 
ment  by  Gutin' s  theory,  there  have  been  a  number  of  outstanding  studies 
(References  1,  7,  14,  17,  30,  34)  involving  propeller  rotational  noise 
generation  and  propagation.  These  studies  have  extended  the  rotational 
noise  theory  developed  by  Gutin  to  account  for  noise  in  the  near  field, 
thickness  noise,  and  noise  field  distortion  due  to  source  translational 
motion.  It  is  desirable  to  make  use  of  the  techniques  developed  through 
these  studies  to  obtain  increased  accuracy  where  Gutin' s  model  is  found 
to  be  inadequate. 

The  study  undertaken  here  has  followed  along  the  lines  of  those 
mentioned  above  in  that  it  seeks  to  extend  Gutin’ s  basic  theory  for  noise 
generated  by  a  particular  type  of  propeller,  in  this  case  the  helicopter 
rotor.  Gutin' s  theory  is  still  commonly  used  for  propeller  noise  predic¬ 
tion  but  is  subject  to  the  following  limitations: 

L.  Uniform  axial  inflow  is  maintained  throughout  the  propeller 
disk. 

2.  Field  points,  at  least  5  diameters  from  the  hub,  are  calculated 
(far  field). 

3.  Only  first  harmonic  noise  is  calculated. 

4.  Propeller  speed  normal  to  the  axis  is  held  below  approximately 
Mach  0. 3. 

In  the  case  of  the  helicopters,  the  inadequacy  of  Gutin' s  theory  is 
obvious;  diameters  are  large  and  inflow  is  nonuniform  because  of  pre¬ 
dominantly  nonaxial  translational  motion.  Most  rotor  noise  field  points 
of  interest  fall  within  the  "near  field",  and  generally  the  first  rotational 
noise  harmonic  falls  below  the  audible  frequency  range,  making  it  less 
important  than  its  harmonics. 

In  general,  the  intensity  of  these  harmonics  cannot  be  predicted 
accurately.  This  report  represents  a  newly  developed  mathematical 
theory  which  removed  the  limitations  of  Gutin' s  theory  as  applied  to  the 
rotational  noise  of  single  rotor  helicopters. 
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The  prediction  method  utilizes  blade  section  loading,  both  steady  and 
varying  at  integral  multiples  of  the  main  rotor  rotational  speed,  as  a 
function  of  radius  and  azimuth.  The  solution  describes  as  many  har¬ 
monics  of  rotational  noise  as  the  input  data  allows  in  both  the  near  and 
far  field.  The  periodically  varying  blade  loading  comprises  the  variable 
inflow  condition.  This  solution  has  been  programmed  by  Sikorsky  Air¬ 
craft  as  an  IBM -7094  procedure. 

DERIVATION  OP  VARIABLE  INFLOW  ROTATIONAL  NOISE  SOLUTION 


The  rotor  system  noise  source  will  be  represented  by  a  surface  of 
stationary  dipole  radiators  which  simulate  the  normal  pressure  distribu¬ 
tion  in  the  rotor  disk.  Shear  forces  as  well  as  thickness  effects  are 
ignored.  It  can  be  shown  (References  7  and  18)  that  for  sound  sources 
in  rectilinear  motion,  the  acoustic  pressure  in  the  sound  field  can  be 
exactly  represented  by  the  sum  of  pressure  from  stationary  sound 
sources  placed  in  the  path  of  the  moving  source  which  radiate  only  when 
the  moving  source  passes  by.  In  this  case  the  stationary  sound  sources 
are  the  dipoles  in  the  rotor  disk  which  radiate  only  when  a  blade 
(corresponding  to  the  moving  source)  passes  by. 

Solution  of  Wave  Equation 


For  the  range  of  rotational  Mach  numbers  and  sound  frequencies 
considered  here,  the  pressure  field  due  to  a  rotating  force  is  the  sum  of 
pressure  fields  from  stationary  dipoles,  and  the  radiation  from  the 
stationary  dipoles  satisfies  the  nonhomogeneous  wave  equation 


2  | 
V  p  — 


=  V-f 7 


Urp 

where  f^is  an  arbitrary  force  per  unit  volume. 


(9) 


Since  this  wave  equation  is  for  radiation  from  a  stationary  force 
into  a  stationary  medium,  it  is  exact  for  only  hover  conditions.  How¬ 
ever,  as  long  as  the  vehicle  translational  Mach  number  is  below  0. 3, 
the  accuracy  of  the  wave  equation  is  considered  adequate,  as  shown  in 
Reference  23. 

Figures  2  and  3  show  the  geometry  of  the  rotor  system  used  in  the 
analysis  and  the  sound  field  coordinate  system.  The  force  vectors 
F(r,V'’#  t)  are  restricted  to  the  region  covered  by  the  rotor  disk,  and  the 
field  point  is  outside  of  a  sphere  having  the  diameter  of  the  rotor  disk 
and  centered  at  the  rotor  hub. 
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The  solution  to  equation  9  for  concentrated  forces  is  given  by 

p(x, 

where 

s  = 

is  the  distance  from  the  elemental  source  point  to  the  field  point. 

Since  the  forces  of  interest  here  are  harmonic  at  integral  multiples 
of  the  rotor  blade  passage,  the  force  vector  may  be  represented  by  a 
Fourier  series: 

P  (r  ,V,t)  *  <r M  e™"u*t  (ii) 

where  the  steady  component  is  neglected  because  it  is  not  an  acoustic 
signal.  The  component  of  the  vector  f  in  the  X  direction  is 

FxO^-f'  3*x(r,v0  eimn*t 

fef  ^ 

where  (fW)  is  the  component, of  (f,f)  in  the  X  direction. 

Similar  components  exist  in  the  y  and*  directions. 

Using  these  expressions  in  equation  10  yields 

a  e‘ mnn|) 

aySL9«ye 

H  D- }  (12) 
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DIRECTION 

OF 

FLIGHT 


Figure  2.  Rotor  System  Geometry  Figure  3.  Field  Point  Location 
which  reduces  to 


This  expression  describes  the  sound  pressure  at  the  field 

point  x,  y,  I  at  timet  due  to  the  force  components'^ 
acting  at  the  point  in  the  rotor  disk.  Now  the  componentsoi  equa¬ 
tion  13  are  defined. 


As  previously  defined 


s  +(y-/ ,/+(*-*■)* 

which  in  spherical  coordinates  (Figure  3)  is 


I  2  2  2 

S  coso"  cos0+rcos?)  +(Rcosrsin0-rsinV^  +  (Rsino) 
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which  reduces  to  _ 

$=^R  +r  -zrRcose-cos(0-l^  (i< 

This  is  an  exact  representaticm  and  holds  for  both  near  and  for  field. 
From  Appendix  111 


s-f  COS  *  CO$ 9 

iy  S 

Equation  15  is  valid  outside  of  one  rotor  diameter  from  the  hub. 


The  differential 


_,„S. 


becomes 


•  •  +  ami 

T  c  J 


-ImnA 


Placing  equations  15  and  16  into  equation  13  yields  for  the  harmonic  Vn 

/  +i  R  -f)r 

n  fy  v  sr  T1  *  ■  11—  g  ro<<r  rncH  4. 


r*  p».*ct,s,rcose  + 

cos<r  sin 6  +  Sin<rj.  +  x| 


szssr 


Blade  Aerodynamic  Loads 


acting  on  the  rotor  blade  remain  to  be 


Blade  pressures  for  this  study  were  obtained  from  flight  test  mea¬ 
surements  taken  on  a  Sikorsky  S-58  (Army  CH-34)  at  NASA,  Langley 
Field,  Virginia  (Reference  29).  This  data  consists  of  differential  pres¬ 
sure  measurements  made  mi  (me  main  rotor  blade  during  various  flight 
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conditions.  Data  was  recorded  at  several  chordal  stations  and  integrated 
over  the  chord  at  15-degree  azimuth  intervals,  yielding  section  loadings 
at  these  points.  This  was  done  at  several  blade  spanwise  stations.  At 
each  station,  a  harmonic  analysis  was  performed  on  the  section  loading- 
azimuth  position  curve  resulting  in  a  steady  plus  10  harmonics  of  the 
loading.  This  harmonic  blade  loading  data  is  utilized  for  the  present 
study. 


Transformation  of  Rotating  Concentrated  Forces  to  Stationary  Forces 


Figure  4.  Blade  Normal  Differential  Pressure  Location 

The  force  acting  cm  a  small  blade  area  is  Apdtdfand  is  considered  to  be 
concentrated  force  when  dTdr  is  very  small.  This  same  force  will  be 
considered  to  be  acting  on  a  differential  area  at  the  rotor  disk  rdrdl^ 
when  the  blade  area  is  covering  the  rotor  disk  area  (Figure  2).  The 
time  history  of  pressure  as  two  blades  pass  over  the  area  rorJf  is 
shown  in  Figure  5. 


Figure  5.  Pressure  Time  History  of  Blade  Passing  Over  an  Elemental 

Disk  Area 
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Note  chat  the  only  time  a  pressure  exists  is  when  a  blade  covers  the  area. 
The  amount  of  time  an  elemental  area  is  covered  is  finite  and  equal  to 
the  time  it  takes  for  the  blade  to  pass  a  given  point.  That  is 
The  time  between  blades  is  equal  to  the  time  required  for  a  blade  to 
traverse  2ir/T\  radians  or  AT*2tr/nA.  The  magnitude  of  die  force  acting 
on  the  rotor  disk  is  aP rArty  and  its  direction  is  normal  to  the  chord. 
After  shifting  scale  and  approximating  the  chordwise  pressure  distribu¬ 
tion  as  a  step  function  (see  Appendix  IV),  the  function  of  Figure  6  is 
obtained. 


Figure  6.  Step  Function 


The  cosine  series  expansion  of  this  function  is 

AP(r.*t)=  f  cos(rtmrjt  -urn  Y^-  _ 


where 


and  is  the  average  chordal  pressure  on  the  blade. 


(18) 

(19) 


Since  average  pressure  data  is  available  (Reference  29)  as  blade 
loading  data  L  (r,Y)  and 

L(r,VJ  =  a  aF 

equation  19  becomes 


*m=  sir  •KW  <»> 
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Substituting  20  into  18  yields: 

Ap(r,r,i)  =  Yjfiri  sinpjgL)L(r,V') co$(mnai-mnr-  (21) 

Putting  this  into  complex  form  and  multiplying  by  irdriV'  to  change  to 
force  yields 


F(r,r,t)=  e<™M-*»+-W)  rJr<lr  (22a) 

w=r 

and  from  Equation  11 

9m  =S§J  MW)  Hr.*)  ;‘nn(r+&)  r  Jr  Jr  (22b) 


Equation  22  corresponds  to  the 
force  exerted  on  the  rotor  disk  by  a 
passing  blade  in  the  direction  nor¬ 
mal  to  the  blade  chord.  From 
Figure  7,  note  that  the  out-of-plane 
component  of  is 

9m  1  =9m  cosP  (23a) 

and  from  Figures  2  and  7  the  in¬ 
plane  components  are 

-3mx~3m  sin0  si"V'  <23b) 
-3my=3msin^ cos*  <23c> 


Figure  7.  Aerodynamic  Forces  where 
Generated  by  a  Blade 

fi  (r~rtt)+Pi  e»s  r+^sin  V 


as  derived  in  Appendix  V. 


Substituting  equation  22  into  equation  23  and  substituting  the  result 


in  equation  17  gives 

pw(xiy,H,t).e 


.imnati 


-^j-|sin/9  smtcos  Ocosr  -sin  /9  cosV'cosrsin$+cosfl  sin*  rdrJfJ- 
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This  expression  after  integrating,  transforming  to  the  uniform  disk  load¬ 
ing  case,  and  introducing  the  small  angle  assumptions  made  in  Refer¬ 
ences  11  and  7,  reduces  to  the  sum  of  equations  21  and  22  of  Reference  7 
for  Mach  zero.  It  represents  the  sound  pressure  at  a  point  x,  y,  Z,  at 
time  t  of  the  jfi  'th  harmonic  of  rotational  noise  due  to  the  blade  section 
loading  L(r,  the  radius  Tof  the  blade  when  the  blade  is  at  azimuth 
location .  The  equation  may  be  integrated  over  the  rotor  disk  to  account 
for  all  sources  and  the  root  mean  square  (rms)  value  of  the  pressure 
with  respect  to  time  may  be  taken  to  correlate  with  sound  pressure 
measurements.  These  operations  result  in  the  following  solution: 


*^sin£  cos«'sin(yr-0)  +  cotfi  sin  irj*  r  dr  kt 
^sinpc«csin(V'-6>)  +cos^  sinrj  r^r 


for  which 


U  a  vnn 


(26b) 


(26c) 


This  integral  solution  (equation  26)  represents  the  rms  sound  pressure 
at  a  point  %,y,M  for  rotational  noise  at  the  m'th  harmonic  of  blade 
passage.  This  solution  may  be  used  for  calculating  the  rms  sound 
pressure  at  any  point  in  the  near  or  for  field  cate  diameter  or  more  from 
the  rotor  hub.  The  number  of  harmonics  which  may  be  calculated  accu¬ 
rately  is  limited  only  by  the  detail  of  the  available  blade  loading  data. 
The  solution’s  accuracy  decreases  with  increasing  rotor  system  trans - 
lational  speed  up  to  Mach  0. 3,  which  is  the  practical  limit  for  accuracy. 


The  solution  has  been  programed  as  an  IBM-7094  procedure  to 
facilitate  calculation.  The  program  is  described  in  Appendix  VI. 
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A  trial  case  was  run  with  steady  loads  only  placed  about  the  0. 8 
radius  point  of  the  rotor  as  assumed  in  Gutin's  derivation  (Reference  11). 


Using  relationships  between  blade  pitch,  loading,  radius,  etc. ,  the 
torque  and  thrust  for  this  case  were  calculated  for  the  S-58  rotor  system. 
The  derivation  is  explained  and  parameters  for  the  example  are  given  in 
Appendix  VII.  Calculations  of  first  harmonic  noise  were  made  at  a  320- 
foot  radius  by  both  the  Sikorsky  program  and  by  Gutin's  equation  for 
various  elevation  angles.  The  results  of  these  calculations  are  tabulated 
in  Table  IV  and  shown  in  Figure  8. 


Figure  8.  Comparison  of  Sikorsky 
Program  and  Gutin  Equation 
Results  at  320- Foot  Radius 


Agreement  between  the  two  solutions 
is  excellent  except  for  the  point 
where  <T  is  -30  degrees.  At  this 
spatial  point,  Gutin's  solution  in¬ 
volves  a  small  difference  between 
two  large  numbers,  and  the  number 
of  significant  digits  used  in  the  hand 
calculation  was  not  sufficient  to 
define  the  true  magnitude. 

To  further  check  the  solution's 
validity,  a  similar  calculation  was 
made  comparing  the  program  output 
with  the  near  field  case  (Figure  4  of 
Garrick  and  Watkins,  Reference  7). 
The  comparison  is  shown  in  Figure 
9  for  zero  Mach  number. 


x/d  IN  FIELD  POINT  PLANE 


Figure  9.  Comparison  of  Results  From  Developed  Solution  and  Garrick 

and  Watkins'  Results 


16 


The  agreement  between  the  calculated  values  and  those  of  Garrick  and 
Watkins  is  excellent,  the  maximum  deviation  being  1  db  up  to  1/2  diameter 
downstream  of  the  propeller. 

EXPERIMENTAL  TEST  PROGRAM 


General 


Measurements  of  rotational  noise  from  an  S' 58  (CH-34)  helicopter 
were  made  at  Bridgeport  Airport  in  Stratford,  Connecticut,  on  April  5, 
1965.  Physical  and  test  parameters  of  the  aircraft  are  listed  in  Table  V. 
Instrumentation  is  described  in  Appendix  VIU.  The  purpose  of  the  test 
was  to  provide  data  for  comparison  with  the  calculated  data  of  Phase  II. 

It  was  found  during  testing  that  although  rotor  noise  harmonics  up  to  the 
tenth  were  sometimes  distinguishable  above  background  noise,  harmonics 
above  the  fourth  (approximately  60  cps)  were  generally  masked  by  engine 
and  tail  rotor  noise. 

Description  of  Test  Helicopter 

The  S-58  (Figure  10)  is  a  single  main  rotor  reciprocating  engine 
helicopter  in  the  13, 000-pound  weight  class.  The  main  rotor  is  4 
bladed,  with  0012  airfoil  shape  and  untapered  square  tip  blades. 


Figure  10.  S-58  Helicopter 
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Acoustical  Measurements 


The  test  layout  is  shown  in  Figure  1. 


Figure  11.  Rotational  Noise  Test  Layout 

The  ambient  wind  was  northeasterly  at  approximately  4  miles  per  hour, 
gusting  occasionally  to  10  miles  per  hour.  The  flight  path  for  cruise 
measurements  was  from  east  to  west  directly  over  the  south  edge  of  the 
runway.  The  ship’s  altitude  was  200  feet  for  all  cruise  and  transient 
measurements. 

For  a  time  and  location  reference,  a  flag  was  dropped  as  the  ship 
passed  the  flag  station,  signalling  the  measurement  station  to  put  a  short- 
duration  400-cycle-per-second  signal  tone  on  the  recording  tape.  Using 
the  aircraft's  indicated  airspeed,  its  position  was  calculated  from  this 
signal.  One-third-octave  band  analysis  was  used  and  the  frequencies 
under  study  went  only  up  to  60  cps,  thereby  eliminating  the  possibility  of 
the  400-cps  tone  interfering  with  the  records. 

Because  the  noise  was  a  series  of  pure  tones,  actual  bandwidth  was 
not  important  in  determining  levels.  One-third-octave  filters  were  used 
so  as  to  define  the  individual  harmonics  and  to  accommodate  the  doppler 
shift  for  the  flyby  conditions. 


18 


Because  of  safety  regulations,  hover  measurements  were  made  at  an 
altitude  of  500  feet.  The  hover  location  was  above  the  south  edge  of  the 
runway. 

Data  was  recorded  with  the  microphone  held  approximately  three 
feet  above  the  ground.  The  runway  between  the  aircraft  and  microphone 
was  hard- surfaced.  In  the  frequency  range  of  interest,  15  to  60  cps,  the 
wavelength  of  75  to  18  feet  precluded  interference  or  reinforcement  be¬ 
tween  the  incident  and  reflected  waves. 

Noise  data  was  recorded  for  hover;  40-knot,  80-knot,  and  110- knot 
cruise;  and  70-knot  cyclic  pullout.  These  conditions  correspond  to  the 
blade  loading  data  of  Reference  28  (Tables  4,  8,  13,  19,  and  111,  res¬ 
pectively). 

CORRELATION  OF  MEASURED  AND  CALCULATED  ROTATIONAL  NOISE 


Correlation  of  measured  and  calculated  rotational  noise  levels  is 
summarized  for  the  first  harmonic  in  Table  VI.  The  calculated  levels 
shown  are  the  maximum  attained  for  flybys  using  NASA  (Reference  29) 
blade  loading  data.  The  measured  levels  are  the  corresponding  maxi- 
mums  for  the  test  flybys.  The  blade  loading  tables  from  Reference  29 
which  were  used  are  listed  in  Table  VII.  The  spread  of  data  for  the 
measured  hover  condition  is  due  to  time  variation.  Levels  calculated  by 
Gutin's  equation,  using  thrusts  and  torques  as  shown  for  the  S-58,  are 
included  for  reference.  Note  that  the  calculated  and  measured  levels 
agree  within  3  db,  while  Gutin's  equation  gives  levels  as  much  as  19  db 
low.  Gutin's  method  is  most  nearly  accurate  (3  db  low)  during  hover 
where  the  inflow  is  almost  uniform.  Obviously,  nonuniform  inflow  has  a 
significant  effect  on  rotor  first  harmonic  noise  levels.  The  effect  of  non 
uniform  inflow  is  even  more  pronounced  for  the  higher  noise  harmonics. 
The  conclusion  is  that  Gutin's  equation  is  not  an  adequate  method  of 
predicting  rotor  noise  levels  during  nonuniform  inflow. 

Figure  23  shows  calculated  and  measured  levels  for  a  500 -foot 
hover.  The  program  calculated  the  measured  data  much  more  closely 
than  Gutin's  equation,  which  predicts  a  rapid  dropoff  in  level  with 
increasing  harmonic. 

It  is  felt  that  the  agreement  between  the  calculated  and  measured 
levels  for  the  higher  harmonics  would  have  been  closer  than  that  shown 
if  the  flight  test  conditions  under  which  the  noise  measurements  were 
made  were  identical  with  those  for  which  the  blade  loading  data  was 
taken.  According  to  Reference  29,  the  blade  loading  test  aircraft  was 
"hovering  in  light  wind",  which  implies  slight,  relatively  steady  wind 
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conditions,  while  during  the  noise  measurement  test,  winds  ranged  from  4 
to  10  knots.  The  lack  of  agreement  between  the  measured  and  calculated 
levels  for  the  higher  harmonics  can  be  attributed  partly  to  the  uncer¬ 
tainty  in  blade  loading  harmonics  caused  by  the  variable  wind  velocity. 

Levels  for  the  cruise  and  transient  flybys  are  shown  in  Figures  24 
through  27.  As  shown,  the  data  is  1-1/2  db  low  at  the  plus  250-foot  and 
minus  250-foot  points.  This  error  results  from  the  frequency  change 
due  to  doppler  effects  reaching  the  band  limits  of  the  one- third-octave 
filter  used  in  the  analysis.  Correlation  of  measured  and  calculated  data 
is  good  for  the  first  two  harmonics.  It  appears  for  some  of  the  cases 
that  the  distance  scale  is  shifted.  This  may  be  attributed  partially  to 
lag  in  the  time  reference  system  used  in  the  test  and  partially  to  uncer¬ 
tainty  as  to  the  test  aircraft's  ground  speed. 

Agreement  between  the  third  and  fourth  harmonics,  and  presumably 
for  all  higher  harmonics,  is  poor.  This  can  be  attributed  to  differences 
in  flight  conditions  between  blade  loading  test  and  noise  test,  as  men¬ 
tioned  earlier,  and  to  the  lack  of  sufficient  higher  harmonic  content  in 
the  NASA  blade  loading  data  (Reference  29).  Blade  loading  containing  a 
greater  proportion  of  higher  harmonics  would  yield  more  nearly  accurate 
higher  harmonics  of  noise  as  explained  in  Phase  III. 

/ 

The  observed  chordal  loading  distribution  for  the  data  used  in  this 
report  conforms  to  the  general  shape  of  Figure  12,  which  yields  a  har¬ 
monic  level  distribution  similar  to  that  of  the  rectangular  distribution 
used.  It  should  be  noted  that  if  the  blade  chordal  pressure  distribution 
had  exhibited  the  form  of  the  (me  shown  in  Figure  13,  the  rectangular 
chordal  pressure  distribution  used  in  the  program  would  yield  much 
lower  noise  harmonic  levels  than  actually  exist.  For  all  common  cases 
of  propeller  or  rotor  loading  where  forces  on  the  disk  act  in  one  direc¬ 
tion,  the  assumption  of  uniform  loading  over  the  chord  should  yield  real¬ 
istic  values  for  calculated  noise.  The  effects  of  variation  in  chordal 
loading  are  treated  in  detail  in  Reference  31. 
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FRACTION  OF  CHORD  FRACTION  OF  CHORD 

Figure  12.  Chordwise  Normal  Force  Figure  13.  Chordwise  Normal  Force 
Distribution  Distribution  With  Reverse  Flow 

APPLICATIONS  OF  THE  ROTATIONAL  NOISE  SOLUTION 


The  theory  derived  here  for  rotational  noise  generation  lends  itself 
to  the  evaluation  of  rotor  design  parameters  better  than  the  classical 
theory  (Reference  11)  for  conditions  of  uniform  inflow.  For  the  more 
significant  condition  of  nonuniform  inflow  it  becomes  a  unique  tool  for 
study.  The  drawback  is  that  the  correct  harmonic  blade  loading  is 
required  as  input,  and  this  is  not  available  for  harmonics  above  the 
fifth.  Analysis  (References  8  and  32)  is  available  to  calculate  these 
blade  loadings;  however,  the  rotor  inflow  (References  21,  24,  and  25)  is 
not  known  in  sufficient  detail  to  allow  the  calculations  to  be  made.  Lack 
of  information  regarding  blade  wake  geometry  is  the  problem  in  predict¬ 
ing  rotor  inflow  characteristics.  The  importance  of  the  need  for  higher 
harmonic  blade  loading  cannot  be  overemphasized,  as  the  usefulness  of 
the  rotational  noise  solution  depends  on  further  progress  in  this  area. 

In  the  case  where  some  blade  loading  harmonics  are  known,  even  if  only 
hand  calculated  steady  loads,  the  distortion  of  the  classical  propeller 
noise  field  {Reference  11)  due  to  planform,  twist,  blade  collective^and 
cyclic  pitch  can  be  calculated  in  both  the  near  and  far  fields.  The  effects 
of  any  other  parameter  changes  are  easily  dealt  with  if  blade  loadings  are 
known.  Note  that  the  solution,  in  its  present  form,  is  not  directly  appli¬ 
cable  to  systems  with  multiple  rotors  in  juxtaposition. 
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PHASE  III 


ROTOR  DESIGN  TECHNIQUES  FOR  ALLEVIATION  OF  NOISE 
VORTEX  NOISE  REDUCTION 

The  research  in  Phase  I  has  shown  how  rotor  vortex  noise  is  gener¬ 
ated  and  propagated.  The  theory  developed  permits  accurate  prediction 
of  the  noise  amplitude  and  frequency  spectrum  directly  from  operational 
and  physical  parameters  of  the  rotor  system.  In  this  section,  the  theory 
and  test  data  of  Phase  I  are  used  to  evaluate  the  amount  of  noise  reduc¬ 
tion  that  can  be  achieved  through  the  proper  choice  of  rotor  system 
design  parameters.  The  formulas  are  valid  for  the  conventional  square- 
tipped  blade  with  uniform  inflow. 

The  work  in  Phase  I  has  pointed  out  the  importance  of  tip  speed  and 
aerodynamic  conditions  at  the  tip  in  noise  generation.  The  intensity  of 
the  tip  vortex  and  unstable  flow  near  the  tip  also  affect  rotor  aerodynamic 
efficiency.  A  company-sponsored  research  program  to  improve  rotor 
efficiency  by  modification  of  the  tip  planform  achieved  increases  in 
efficiency  and  reductions  in  noise.  Reduction  of  tip  vortex  induced 
velocity,  by  approximating  an  elliptical  load  distribution, was  effective  in 
reducing  noise. 

Effect  of  Blade  Area,  Tip  Speed  and  Lift  Coefficient 

The  theoretical  and  experimental  results  of  Phase  I  have  shown  that 
the  most  important  parameters  governing  noise  generation  are  blade 
area,  tip  speed  and  lift  coefficient  (or  angle  of  attack).  The  relationship 
can  be  formalized  as  follows: 

I“j^  K(Vr)*AB(CL)2  (27) 

The  parameters  Vp,  C|_,  and  Ag  also  govern  the  thrust  of  the  rotor 
system,  as  can  be  seen  from  the  following  basic  relationship: 

T=±rV2cLAB  (28) 

Although  the  above  equation  assumes  a  two-dimensional  section  with 
velocity  V  for  the  purpose  of  noise  estimation,  we  can  let  the  velocity  V 
be  that  at  the  0. 7  radius  and  assume  that  Ci  is  the  lift  coefficient  at  a 
mean  angle  of  attack.  u 

Since  the  same  three  basic  parameters  appear  in  the  equations  for 
noise  generation  and  lift,  trade-offs  between  these  parameters  for  noise 
control  will  directly  affect  lift.  However,  the  appearance  of  blade 
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velocity,  V,  to  the  sixth  power  in  the  noise  equation  allows  for  limited 
reduction  of  tip  speed  without  sacrificing  lift. 

In  Figure  28,  calculated  overall  vortex  noise  levels  are  shown  for  a 
range  of  tip  speeds  with  corresponding  blade  areas  ;.o  maintain  lilt.  The 
curves  shown  are  for  coefficients  of  lift  of  0. 4  and  0. 2.  These  curves 
are  for  thrust  of  20, 000  pounds.  The  effect  cm  vehicle  weight  of  the 
changes  in  blade  areas  necessary  for  practical  rotor  system  designs 
was  not  considered. 

The  curves  show  that  sharp  decreases  in  overall  vortex  noise  can  be 
achieved  by  reduction  of  the  tip  Mach  number.  However,  as  mentioned 
above,  the  larger  and  heavier  blade  systems  have  to  produce  more 
thrust  to  maintain  the  same  payload.  The  noise  associated  with  the 
increased  thrust  partially  offsets  the  reduction  shown  in  Figure  28. 

The  effect  of  thrust  on  the  overall  vortex  noise  is  shown  in  Figures 
29  and  30.  The  slope  of  the  curves  is  approximately  20  times  the  log  of 
the  thrust  ratio.  Because  of  the  wider  variation  in  tip  speeds,  the 
CH-53A  system  curves  show  more  clearly  that  tip  speed  is  more  impor¬ 
tant  than  angle  of  attack  in  reducing  the  vortex  noise.  For  equal  thrusts 
and  blade  area,  reducing  the  tip  speed  necessitates  increasing  the  angle 
of  attack.  However,  the  reduction  in  noise  level  due  to  lowering  tip 
speed  from  800  to  630  feet  per  second  is  about  5  db,  even  with  this 
compensating  increase  in  angle  of  attack. 

Similar  curves  can  be  constructed  for  any  desired  thrust  by  using 
the  vortex  noise  equation  developed  in  Phase  1  (equation  6).  Selection  of 
rotor  design  parameters  can  be  made  with  respect  to  noise  and  reconciled 
with  other  mission  requirements. 

Elliptical  Tip  Loading 

Although  a  limited  amount  of  vortex  noise  reduction  can  be  achieved 
by  trade-offs  between  tip  speed,  blade  area  and  thrust  coefficient,  a 
more  fruitful  area  of  noise  reduction  was  discovered  during  tests  of 
blades  having  elliptical  blade  loading  at  the  tip.  Figure  31  shows  the 
planform  for  a  trapezoidal  tip  and  the  standard  square- tipped  blades 
studied  in  the  remainder  of  this  report.  The  trapezoidal  shape  was 
adopted  as  an  approximation  of  the  elliptical  blade  loading.  It  was 
expected  that  elliptical  loading  would  reduce  the  induced  velocity  of  the 
tip  vortex  and  thereby  increase  lifting  efficiency. 

At  the  onset  of  the  program  it  was  anticipated  that  some  noise  reduc¬ 
tion  would  result.  The  noise  reduction  proved  to  be  substantial  and 
appeared  to  be  directly  related  to  increased  efficiency. 
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In  Figure  32,  the  overall  vortex  level  is  plotted  against  load  for  three 
tip  speeds.  The  tip  speeds  are  the  same  as  those  in  Figure  29.  Except 
for  the  tip  change,  the  rotor  system  is  identical  to  that  plotted  in  Figure 
29.  The  reference  line  from  Figure  29  is  repeated  to  allow  direct  com¬ 
parison  of  the  levels.  It  is  interesting  to  note  that  the  greatest  difference 
between  the  square  and  elliptical  tip  noise  levels  is  at  low  thrust.  The 
sound  pressure  levels  of  the  elliptical  tips  increase  more  rapidly  with 
increased  thrust.  At  high  thrust  the  difference  between  the  two  designs 
is  much  smaller.  At  the  normal  disk  loading  of  about  19, 000  pounds  for 
this  size  rotor,  the  reduction  in  overall  noise  is  about  7  db. 

Figure  33  shows  the  spectra  for  both  square  and  trapezoidal  tip 
blades  at  the  extremes  of  the  usual  range  of  angles  of  attack.  The 
spectra  for  the  trapezoidal  tips  show  much  less  noise  in  the  octave  con¬ 
taining  the  Strouhal  frequencies.  Each  of  the  elliptical  spectra  showed  a 
rise  in  the  last  octave  which  is  unexplained.  The  square -tipped  blade 
in  Figure  33  is  evidently  experiencing  stall,  as  discussed  in  Phase  I, 
because  of  the  secondary  peak  in  the  1200  to  2400  octave.  For  the  same 
pitch  angle  the  trapezoidal  tip  stills  appears  to  be  below  the  onset  of 
stall. 

Figure  34  shows  the  relative  levels  of  the  300  to  600  octave  which 
contains  the  Strouhal  frequencies  over  a  range  of  loads.  The  reduction 
in  the  vortex  noise  at  the  Strouhal  frequencies  is  obvious.  Also  it  can 
be  seen  that  above  19, 000  pounds  the  noise  in  this  key  octave  increases 
more  rapidly  with  load  for  the  trapezoidal  tip. 

The  reduction  in  tip  vortex  strength  deduced  from  the  reduction  of 
Strouhal  frequency  noise  is  also  demonstrated  by  less  modulation  of  the 
noise  amplitude  as  shown  in  Table  VIII.  Modulation,  which  results  from 
movement  of  the  source  relative  to  the  point  of  observation,  depends  on 
the  strength  of  the  moving  source.  If  the  moving  source  is  much 
stronger  than  the  average  over  the  field,  the  modulation  is  greater. 

Recommendations 


Main  rotor  vortex  noise  can  be  reduced  by  reducing  rotor  tip  speed 
and  by  improving  load  distribution  at  the  tip.  It  is  recommended  that 
minimum  tip  speed  consistent  with  good  aerodynamic  design  be  used. 

In  addition,  reduction  of  the  induced  velocity  of  the  tip  vortex  (core 
thickening)  will  result  in  substantial  vortex  noise  reduction  without  loss 
in  efficiency.  Vortex  noise  reductions  of  from  7  to  10  db  have  been 
achieved  with  this  approach.  Since  vortex  noise  falls  in  a  frequency 
range  for  which  hearing  acuity  is  high,  it  can  be  important  in  determin¬ 
ing  detectability,  as  can  be  seen  in  Figure  35.  The  tactical  potential  of 
good  blade  tip  design  is  therefore  apparent. 
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ROTATIONAL  NOISE  REDUCTION 


Reduction  of  Noise  Due  to  Steady  Blade  Loadings 

Far  field  rotational  noise  due  to  steady  blade  loading  and  its  relation¬ 
ship  to  the  physical  parameters  of  propeller  systems  are  described 
mathematically  by  Gutin  (Reference  II).  This  relationship  has  been  put 
into  graphic  and  tabular  form  (References  3,  15, and  16)  which,  in 
combination  with  aerodynamic  and  weight  considerations,  makes  the 
matter  of  parameter  selection  for  low  noise  straightforward.  Additional 
investigations  (References  12,  l$>and  27)  have  defined  the  trade-offs 
necessary  to  obtain  lower  noise  levels.  These  studies,  in  many  cases, 
have  concentrated  on  the  propeller  noise  as  it  affects  aircraft  detect¬ 
ability.  Aircraft  detectability  relates  noise  generated  by  an  aircraft  to 
the  ability  of  a  human  to  detect  its  presence  audibly.  Detectability  is 
extremely  important  in  military  applications,  as  an  aircraft's  acoustic 
detectability  is  a  factor  in  determining  combat  effectiveness  and  proba¬ 
bility  of  survival.  The  concept  of  detectability  will  be  used  as  a  basis 
for  evaluation  of  rotor  noise  alleviation  techniques. 

Figure  35  from  Reference  20  shows  the  minimum  detectable  level 
of  aircraft  noise  in  a  low  noise  background.  The  level  up  to  approxi¬ 
mately  250  cps  is  the  threshold  of  hearing  (Reference  6).  Compare  this 
with  an  S-58  helicopter  noise  spectrum  shown  in  Figure  36.  It  is 
obvious  mat  the  fundamental  (normally  between  15  and  20  cps)  of  the 
main  rotor  blade  passage  noise  would  have  to  be  at  a  much  higher  level 
than  the  higher  harmonics  in  order  to  determine  detectability.  In 
general,  the  fundamental  component  of  rotor  rotational  noise  is  below 
audible  level  and  is  not  heard.  Instead,  blade  passage  harmonics  and 
modulated  vortex  frequency  noise  are  heard  and  incorrectly  identified  as 
blade  passage  frequency  noise.  It  is  recognized  then  that  die  fundamental 
of  rotational  noise  of  helicopter  rotor  systems  is  not  important  in  deter¬ 
mining  detectability.  While  Gutin  theory  does  not  accurately  predict  the 
level  of  the  harmonics,  it  is  adequate  for  estimating  the  relative  effect 
of  design  parameters. 

Reduction  of  Noise  from  Nonuniform  Rotor  Loading 

This  study  is  based  on  the  assumption  that  first  harmonic  rotational 
noise  is  not  significant  in  judging  loudness,  and  that  harmonic  levels 
resulting  from  steady  inflow  to  the  rotor  will  vary  directly  as  the  funda¬ 
mental  for  which  noise  control  methods  are  already  known.  To  mini¬ 
mize  the  acoustic  signature  of  a  helicopter,  techniques  must  be  developed 
to  reduce  the  harmonic  content  cf  rotational  noise. 
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According  to  Gutin,  the  noise  level  of  the  individual  harmonics 
decreases  rapidly  with  harmonic  number,  as  shown  in  curve  A  of  Figure 
37.  Gutin' 8  equation  assumes  that  all  blade  loading  is  concentrated  at  a 
single  point  on  the  blade.  This  loading  is  simulated  in  curve  A.  When 
the  concentrated  loading  is  spread  out  over  a  larger  segment  of  the 
blade,  the  upper  harmonics  increase  in  level.  This  implies  that  concen¬ 
tration  of  blade  loading  over  a  smaller  disk  area,  which  might  be  accom¬ 
plished  by  reducing  the  diameter,  for  example,  results  in  lower  harmonic 
noise  levels.  This  relationship  cannot  be  utilized,  however,  because  the 
harmonic  dropoff  gained  by  concentrating  the  load  would  be  more  than 
offset  by  the  increase  in  overall  noise  level  from  the  higher  blade  loading. 

The  fact  that  redistribution  of  steady  blade  loadings  cannot  be  uti¬ 
lized  to  reduce  harmonic  noise  levels  leaves  nonuniform  loading  as  the 
controlling  factor  in  determining  harmonic  level. 

The  influence  of  steady  plus  4-per-rev  blade  loading  on  harmonic 
noise  levels  is  shown  in  Figure  38.  Three  different  levels,  of  4-per-rev 
plus  steady  loading  were  used  to  calculate  the  harmonic  noise  levels  for 
the  S-58  rotor  system.  The  levels  were  calculated  by  the  methods  of 
Phase  II.  The  rotor  system  parameters  and  loads  for  this  hypothetical 
system  are  presented  in  Appendix  IX.  The  harmonic  dropoff  is  substan¬ 
tially  affected  by  the  addition  of  the  variable  blade  loading.  For  example, 
the  second  noise  harmonic  is  now  higher  than  the  fundamental  and  all  the 
levels  of  the  harmonics  are  increased  over  those  generated  by  the  steady 
loading. 

The  4-per-rev  loading  contributes  more  to  the  level  of  the  higher 
harmonics  than  the  fundamental  does.  Calculations  were  made  to  deter¬ 
mine  whether  a  pattern  exists  which  describes  relationship  between 
blade  loading  and  rotational  noise  levels.  The  results  are  shown  in 
Figures  39  to  45,  which  represent  the  same  information  as  Figure  38 
except  that  the  noise  due  to  steady  loading  has  been  removed  to  show 
only  the  result  of  harmonic  loading.  The  increase  in  sound  pressure 
level  caused  by  doubling  the  harmonic  loading  (from  1  /4  to  1/2  of  the 
steady  loading  amplitude)  is  not  always  the  6  db  which  would  result  from 
interaction  between  the  harmonic  and  the  steady  loadings.  The  shapes 
of  the  curves  change  considerably  between  one  blade  loading  harmonic 
and  the  next.  The  higher  blade  loading  harmonics  have  a  greater  effect 
chi  the  higher  noise  harmonics  and  less  effect  on  the  lower  noise  har¬ 
monics.  Notice  that  as  the  blade  loading  harmonic  increases,  the  noise 
produced  by  a  given  level  of  loading  greatly  increases.  In  other  words, 
the  efficiency  of  conversion  from  blade  loading  to  noise  increases  with 
blade  loading  harmonic.  This  effect  is  shown  for  blade  loadings  up  to 
10 -per -rev  in  Figures  46  and  47,  where  noise  level  is  plotted  against 
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rotational  noise  harmonic  for  several  harmonic  blade  loading  conditions. 
This  is  done  for  a  harmonic  loading  level  of  one -quarter  the  amplitude 
of  the  steady  component  plus  the  steady  loading,  but  with  the  noise  attri¬ 
buted  to  the  steady  load  omitted. 

To  get  a  better  idea  of  how  the  noise  producing  efficiency  of  blade 
loading  varies  with  noise  and  loading,  the  data  of  Figures  46  and  47  is 
cross  plotted  in  Figure  48.  The  fact  that  certain  blade  loading  harmonics 
generate  noise  most  efficiently  at  certain  noise  harmonic  frequencies  is 
clearly  shown  here.  Again  the  harmonic  blade  loading  level  used  here  is 
1  /4  that  of  the  steady  loading  and  the  noise  attributed  to  the  steady  load¬ 
ing  itself  has  been  removed.  Even  when  the  noise  levels  due  to  the 
steady  loading  are  added  to  those  of  Figure  48,  as  shown  in  Figure  49, 
it  is  quite  easy  to  produce  a  spectrum  where  many  harmonic  levels  are 
higher  than  the  fundamental.  The  acoustical  efficiency  of  the  many 
blade  loading  harmonics  is  the  cause  of  the  high  levels  of  rotor  and 
propeller  noise  harmonics. 

If  a  helicopter  is  being  designed  for  low  detectability,  it  is  essential 
to  control  the  harmonics  of  blade  loading.  The  more  uniform  the  inflow, 
the  lower  are  the  blade  loading  harmonics.  Although  rotor  system  dyna¬ 
mics  and  aeroelastic  characteristics  play  a  role  in  determining  the  in¬ 
flow  pattern,  the  most  important  acoustic  consideration  is  the  presence 
or  absence  of  wake  interaction  from  other  blades  or  rotor  systems. 

Quantitative  analysis  of  the  pressure  field  generated  by  a  rotor  and 
its  blade  loading  is  difficult.  Improved  correlation  between  measured 
and  predicted  blade  loading  has  been  obtained  by  replacing  steady  'aero¬ 
dynamic  inflow  theory  with  variable  induced  inflow  theory  (References  8 
and  31).  Further  work  is  required,  however,  to  generate  harmonics  of 
higher  order  analytically  with  any  degree  of  accuracy.  Figures  39  to 
45  show  the  importance  of  the  higher  order  blade  loading  harmonics  in 
determining  harmonics  of  rotational  noise.  As  an  example.  Figure  46 
shows  that  if  a  2-per-rev  blade  loading,  of  1/4  the  amplitude  of  the 
steady,  were  present  in  addition  to  the  steady,  the  fourth  harmonic  of 
rotational  noise  would  be  increased  by  7  db.  If  there  were  a  5-per-rev 
blade  loading  of  1/4  the  magnitude  of  the  steady  in  addition  to  the  2-per- 
rev  and  steady  loadings,  the  level  of  the  fourth  harmonic  of  rotational 
noise  would  be  further  increased  by  41  db.  For  a  rotor  with  n  number  of 
blades,  it  appears  as  though  the  with  harmonic  of  rotor  noise  is  most 
sensitive  to  the  Wl(h-l)  loading  harmonic.  To  determine  higher  order 
noise  harmonics  accurately,  it  appears  as  though  there  is  a  need  for 
loading  data  of  much  higher  harmonic  content  than  presently  is  available 
from  measurements  or  analysis. 
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Design  Recommendation 

It  is  concluded  from  the  Phase  IQ  study  that  it  is  necessary  to  con¬ 
trol  rotor  harmonic  loading  in  order  to  control  rotational  noise.  Means 
are  not  presently  available,  however,  to  predict  the  nonuniform  portion 
of  rotor  loading  to  the  accuracy  required  for  use  in  the  rotor  noise 
program.  Gutin's  theory  of  Reference  11  is  adequate  to  predict  trends 
only  for  first  harmonic  rotational  noise  under  steady  inflow  conditions. 
The  charts  of  References  3,  15,  and  16  then  become  useful  in  determi¬ 
ning  the  acoustic  effect  of  blade  parameter  changes,  but  not  die  absolute 
level.  For  predicting  rotational  noise,  for  even  relatively  uniform  in¬ 
flow  conditions,  die  analysis  presented  in  Phase  II  should  be  used.  The 
required  input  data  for  such  predictions  must  include  up  to  at  least  the 
Wl  times  n  harmonic  of  blade  loading.  Here  II  is  the  number  of  blades 
and  HI  is  the  desired  harmonic  of  rotational  noise.  Since  die  state  of  the 
art  does  not  allow  prediction  of  higher  than  the  fourth  or  fifth  harmonic 
of  blade  loading  little  can  be  concluded  as  to  the  parametric  effects  of 
design  changes  cm  harmonics  of  rotational  noise. 

In  addition  to  the  general  requirements  of  low  tip  speed  and  blade 
loading,  it  is  recommended  that  the  blade  tip  be  designed  so  that  loading 
is  distributed  as  evenly  as  possible  and  local  turbulence  is  held  to  a 
minimum.  This  is  the  same  requirement  as  that  specified  for  vortex 
noise  control  and  can  be  accomplished  with  tips  such  as  the  trapezoidal 
variety. 

BLADE  SLAP 


Discussion 


Blade  slap  is  the  sharp  popping  sound  produced  by  a  helicopter  rotor 
during  certain  flight  conditions.  This  characteristic  was  extensively 
investigated  in  this  study  to  define  the  mechanism  of  blade  slap  and 
means  of  controlling  it. 

Results  of  this  study  and  that  described  in  Reference  2  indicated 
that  blade  slap  consists  of  high-amplitude  rotational  noise  plus  highly 
modulated  vortex  frequency  noise.  The  slapping  or  popping  noise  occurs 
when  the  blade  has  rotated  approximately  270  degrees  from  the  tail 
during  forward  flight.  Reference  5  indicates  that  blade  slap  is  due  to 
amplitude  modulation  of  broad-band  noise  during  stall.  It  can  be  seen 
in  Figure  33  and  in  Figure  2  of  Reference  33  that  vortex  noise  levels  are 
a  minimum  of  7  db  higher  for  the  stalled  (separated  flow)  condition  than 
for  the  unstalled  condition. 
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The  contention  that  blade  slap 
is  associated  with  blade  stall  can  be 
demonstrated  as  follows.  Figure  14 
shows  the  angles  of  attack  around 
the  disk  of  a  single  rotor  helicopter 
during  100 -knot  forward  flight.  The 
maximum  angle  of  attack  occurs  at 
an  azimuth  angle  between  270  degrees 
and  360  degrees.  As  speed  is  in¬ 
creased,  the  angle  of  attack  in  this 
area  is  increased.  Ultimately,  stall 
will  occur  at  the  tip  and  will  spread 
toward  the  rotor  hub.  The  occur¬ 
rences  of  stall  and  angle  of  attack 
change  are  explained  quite  clearly  in 
Reference  9.  These  variations  in 
angle  of  attack  are  typical  for  a 
Figure  14.  Single  Rotor  Helicopter  single  rotor  helicopter.  Experi- 
Typical  Angle  of  Attack  Distribu-  mental  data  from  Reference  10  on  a 
tion  rotor  system  experiencing  stall 

(Figure  15)  shows  that  calculated 

results  (dotted  line)  are  reasonably  accurate  and  that  the  region  around 
270  degrees  is  indeed  the  area  which  is  most  susceptible  to  stall  during 
high  speed  flight.  Since  both  blade  slap  and  stall  occur  near  the  270 
degrees  azimuth  position,  it  can  be  concluded  that  slap  is  related  to 
blade  stall  for  this  flight  condition. 

Pilot  reports  state  that  the  blade  slap  occurs  with  high  blade  loading 
such  as  that  encountered  in  high  speed  forward  flight  or  in  a  heavily 
laden  condition.  They  also  report  slap  when  making  powered  descent 
and  in  transition  to  autorotation.  Where  stall  is  imminent,  a  perturba¬ 
tion  in  the  aerodynamic  inflow,  such  as  the  downwash  of  another  rotor 
or  an  encounter  with  a  shed  vortex  from  a  preceding  blade,  may  be 
sufficient  to  induce  the  stall. 
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FORWARD 


Figure  15.  Occurrence  of  Stall  in  the  Rotor  Disk 


With  tandem  rotor  helicopters,  the  lower  rotor  encounters  the 
downward  induced  flow  along  with  shed  vortices  in  a  portion  of  its  disk 
area.  The  intersection  of  these  vortices  and  the  resulting  sudden 
change  in  angle  of  attack  and  flow  separation  cause  nearly  continuous 
slapping.  A  comparison  of  the  levels  of  a  single  rotor  (S-61)  helicopter 
and  a  tandem  rotor  (V-107)  helicopter  during  approach  to  landing  is 
shown  in  Figure  16.  The  ships  are  of  approximately  18, 000 -pound 
design  gross  weight  and  are  performing  the  same  maneuver.  During  the 
maneuver,  both  aircraft  are  producing  impulse  type  noise. 
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OCTAVE  BAND  CENTER  FREQUENCIES 
IN  CYCLES  PER  SECOND 
DISTANCE  70  FEET  FROM  BOTH  SHIPS 

Figure  16.  Single  versus  Tandem  Rotor  Noise  in  Approach 

(Reference  34) 

The  frequency  content  of  the  noise  for  the  two  helicopters  is  similar. 
The  considerably  higher  noise  in  the  20 to  75-cps  octave  indicates  that 
the  tandem  helicopter  experiences  corresponding  changes  in  lift  in  the 
overlap  region.  As  was  shown  in  Phase  II,  the  periodic  change  in  lift  or 
blade  pressure  causes  high  level  noise  at  the  blade  passage  frequency 
and  its  harmonics.  Although  momentary  tip  stall  may  occur  on  the 
single  rotor  ship,  the  fluctuation  in  lift  is  less  intense.  Consequently, 
it  generates  little  more  noise  than  in  the  unstalled  condition. 

Tests  of  blades  with  square  and  trapezoidal  tips  were  made  by 
Sikorsky  Aircraft.  The  trapezoidal  tips  are  designed  to  reduce  the 
induced  velocity  of  the  tip  vortex  from  that  of  the  square-tipped  blade. 
This  lessens  the  depth  of  vortex  noise  modulation  and  also  leaves  less 
of  a  disturbance  in  the  flow  field  to  interact  with  the  other  blades.  The 
differences  between  noise  spectra  of  the  two  blade  types  are  shown  for 
various  power  and  speed  conditions  in  Figures  33  and  34.  Details  of 
the  test  are  explained  in  Phase  I. 
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The  phenomenon  of  compressible  drag  divergence,  which  occurs  with 
stall,  is  also  related  to  blade  slap.  When  drag  divergence  occurs,  drag 
increases  and  power  is  consumed.  The  increase  in  drag  is  due  to  the 
unstable  formation  of  local  shock  waves  which  transform  energy  into 
heat  and  pressure  pulses.  If  sufficiently  strong,  these  shock  waves 
contribute  to  the  blade  slap  noise. 


The  Phase  II  calculation  procedure  is  completely  general,  and  as 
such  will  predict  the  harm  exile  content  of  any  condition,  including  blade 
slap.  Blade  slap  was  not  detected  during  the  measurement  program. 

In  order  to  determine  whether  blade  slap  could  be  predicted  by  the  Phase 
II  analysis,  it  was,  therefore,  necessary  to  simulate  a  blade  loading 
condition  which  would  induce  blade  slap.  Since  blade  slap  occurs  during 
blade  stall,  blade  pressure  data  recorded  during  the  occurrence  of 
blade  stall  was  used  as  a  basis  for  the  harmonic  loading  input.  Blade 
pressure  time  histories  at  85  percent  span  for  the  80- knot  l/2g  turn  are 
shown  in  Figure  17.  From  analysis  of  the  pressure  time  history  at  the 
trailing  edge,  it  was  evident  that  stall  was  present.  This  trace  was 

analyzed  for  harmonic  content  over 
one  blade  revolution.  The  resulting 
harmonics  were  used  as  blade  load¬ 


ing  for  the  entire  chord  at  the  85 
percent  span  of  the  rotor  system  and 
the  noise  harmonics  were  calculated. 
Noise  levels  were  also  calculated 
using  only  the  steady  component  of 
the  blade  loading.  The  resulting 
harmonic  distributions  are  shown  in 
Figure  18  along  with  the  harmonic 
distribution  resulting  from  steady 
80-knot  flight  on  the  same  S-58.  The 
steady  loading  spectrum  exhibits  the 
usual  rapid  dropoff  in  harmonic 
noise  level  associated  with  steady 
loads  alone.  The  80-knot  cruise 


izMmiAMu -m  spectrum  levels  drop  off  also, 

although  not  so  rapidly.  Finally,  the 
Figure  17.  S-58  Blade  Differential  spectrum  generated  using  the  impul- 
Pressure  Time  History  During  An  sive  airloads  experienced  during 
80-Knot,  1/2  g  Turn  at  85  Percent  blade  stall  shows  extremely  high 
Span  -  A  Condition  harmonics. 


of  Partial  Blade  Stall 
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SOUND  PRESSURE  LEVEL 


123  456789  10 


ROTOR  NOISE  HARMONIC 

Figure  18.  Rotational  Noise  Harmonic  Content 

This  spectrum  is  similar  in  shape  to  the  blade  slap  spectrum  measured 
on  an  SH-3A  helicopter  during  an  abnormally  rough  flare  from  level 
flight  to  hover,  at  which  time  the  blades  were  experiencing  localized 
blade  stall.  Although  the  spectrum  shown  is  limited  to  150  cps,  blade 
slap  extends  over  a  frequency  range  of  20  to  over  1000  cps,  and 
increases  vortex  noise  levels  as  well  as  the  rotational  noise  harmonic 
levels.  Reference  1  also  substantiates  this  point.  Thus,  if  detailed 
blade  surface  pressures  are  available,  the  blade  slap  rotational  har¬ 
monic  spectrum  can  be  predicted.  As  shown  in  Figure  17,  stall  can 
be  a  localized  phenomenon,  requiring  detailed  blade  surface  pressures 
to  be  described  adequately.  These  prediction  methods,  however,  are 
not  presently  available.  It  appears  that  the  best  basis  upon  which  to 
judge  the  likelihood  of  blade  slap  occurrence  is  aerodynamic  stall 
criteria. 

Design  Recommendations 

For  minimal  rotational  noise  and  blade  slap,  the  obvious  design 
guidelines  are  to  use  the  lowest  possible  tip  speed  and  to  use  shaped, 
twisted  blades  which  distribute  tip  loading  and  substitute  sheet  vortices 
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I  for  the  stronger  tip  vortices.  Although  blade  stall  appears  to  be  the 

|  major  cause  of  blade  slap,  compressible  drag  divergence,  i.  e. ,  Mach 

I  effects,  cannot  be  ignored.  When  blade  elements  are  in  stall,  com- 

:  pressible  drag  divergence  is  also  present.  To  make  blade  design 

f  recommendations  for  reduced  blade  slap  strictly  on  the  basis  of  reduced 

stall  would  be  premature  at  this  time,  since  blade  design  for  reduced 
tendency  to  stall  may  well  increase  the  extent  of  drag  divergence. 
Further  study  is  therefore  needed  to  define  the  relationships  between 
blade  slap  and  the  degree  of  stall,  drag  divergence,  and  tip  vortex 
induced  velocity  if  positive  blade  design  recommendations  are  to  be 
made. 
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TABLE  I 


CALCULATED  VORTEX  NOISE  FREQUENCIES  AND  THEIR  RELATIVE 
INTENSITY  FOR  THE  CH-53A  ROTOR  BLADE  AT  30, 000  POUNDS  THRUST 


r/R  <x 

(degrees) 

M 

V 

(ft.  sec) 

K 

(ft) 

f* 

(cps) 

Down  Fran 
Tip  Level 
(db-db  at  Tip) 

0.300 

3.31 

0. 187 

209 

0.367 

154 

-31.3 

0. 450 

3.54 

0. 280 

313 

0. 375 

225 

-20.9 

0.600 

3.44 

0. 374 

418 

0.371 

300 

-13.3 

0.700 

3.30 

0. 436 

487 

0.366 

360 

-  9.3 

0.750 

3.20 

0.467 

522 

0. 362 

390 

-  7.5 

0.800 

3. 10 

0. 499 

556 

0.359 

420 

-  5.9 

0.850 

3.00 

0. 530 

591 

0. 355 

450 

-  4.3 

0.900 

2.84 

0. 561 

626 

0.349 

480 

-  2.8 

0.950 

2.64 

0. 592 

661 

0.341 

520 

-  1.3 

0.975 

2. 45 

0.608 

679 

0. 334 

550 

-  0.7 

1.000 

0.00 

0.624 

696 

0.000 

- 

0.0 

.  f.  StV 

w 

where  S-f 

=  0.28 

TABLE  II 


OVERALL  VORTEX  NOISE  LEVELS 


Rotor 

Speed 

<rpm) 

Thrust 

(pounds) 

SPL  Calculated 
by  Equation  6 
<db*) 

SPL 

Measured 

<db*) 

SPL  Calculated 
by  Equation  5 
<db*> 

CH-3C  ROTOR  SYSTEM 

183.0 

13,400 

78 

76 

79 

16,200 

80 

79 

18,700 

81 

81 

20,500 

82 

82 

203.0 

16,300 

82 

81 

82 

18, 100 

83 

•  81 

19,900 

83 

82 

21,400 

83 

83 

213.0 

14,500 

81 

81 

83 

18,200 

83 

83 

20,000 

83 

83 

CH-53A  ROTOR  SYSTEM 

166.0 

24,000 

82 

80 

83 

28,400 

83 

81 

32,000 

84 

83 

36,200 

85 

83 

39,000 

86 

85 

185.5 

25,000 

83 

83 

85 

30, 100 

85 

83 

36,200 

86 

85 

41,600 

87 

87 

215.0 

23, 700 

84 

85 

89 

29,600 

86 

86 

37,900 

88 

89 

43,520 

89 

90 

*  Ref  =  0. 0002  dyne/sq.  cm. 


36 


TABLE  III 


PHYSICAL  CHARACTERISTICS  OF  TEST  ROTOR  SYSTEMS 


CH-3C 

CH-53A 

Number  of  blades 

5 

6 

Diameter  (feet) 

62 

72 

Blade  area  (ft^) 

217 

368 

Solidity 

0.078 

0.115 

Disk  area  (ft^) 

3020 

4070 

Airfoil 

0012 

0012  Modified 

Chord  (feet) 

1.50 

2.16 

Rotor  speed  (RPM) 

183 

166.0 

203 

185.5 

213 

215.0 

Tip  speed  (ft/sec) 

595 

625 

661 

696 

692 

810 

Tip  Mach  number 

0.  532 

0.591 

0.586 

0.625 

0.613 

0.726 
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TABLE  IV 

COMPARISON  OF  FIRST  HARMONIC  NOISE 
CALCULATIONS  BY  GUTIN’S  EQUATION  AND  BY  THE 

SIKORSKY  PROGRAM 


<r* 

degrees) 

Sikorsky 

Gutin 

Sound  Pressure 
(psi) 

SPL 

<db**) 

Sound  Pressure 
(psi) 

SPL 

(db**) 

-45 

1.6x  10'5 

74.8 

1. 58  x  10~5 

74.7 

-30 

2. 92  x  10"5 

80.0 

2. 82  x  10-5 

79.8 

-15 

3.  41  x  10"5 

81.4 

3. 29  x  10“5 

81.0 

0 

2.  69  x  10' 5 

79.3 

2. 58  x  10“5 

79.0 

15 

1. 32  x  10  5 

73.2 

1. 24  x  10’6 

72.6 

30 

3.  39  x  10'6 

61.4 

1. 87  x  10" 6 

56.2 

*  Elevation  angle  =  0°  in  rotor  plane,  positive  above,  negative  below 
**  Ref  =  0. 0002  dyne/ sq.  cm. 


TABLE  V 

FLIGHT  TEST  PARAMETERS 


Test  weight- pounds  approximately  12, 000 


Number  of  blades  4 

Blade  radius -feet  28 

Blade  twist-degrees  -8 

Blade  chord -feet  1. 367 

Test  rotor  speed-rpm  212 

Test  engine  speed-rpm  2400 

Test  rotor-angular  velocity-radians/second  22. 2 
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TABL  E  VI 

MEASURED  AND  CALCULATED  FIRST  HARMONIC  NOISE  LEVELS 
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TABLE  VII 


BLADE  LOADING  TABLES  USED  IN 
PROGRAM 


Blade  Loading  Data  Table  * 

Noise  Measurement  Condition 

4 

Hover  out  of  ground  effect 

8 

40-Knot  cruise 

13 

80-Knot  cruise 

19 

110-Knot  cruise 

111 

70-Knot  cyclic  pullout 

*  Reference  28 
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TABLE  VIII 


17-CY CLE-PER-SECOND  AMPLITUDE  MODULATION:  COMPARISON  OF 
CH-3C  TRAPEZOIDAL  TIP  VERSUS  SQUARE  TIP  BLADES 


Octave 

Average  Depth  of  Modulation  -  db* 

Square  Tip  Trapezoidal  Tip  Difference 

150-300 

8 

5 

3 

300-600 

10 

6 

4 

600-1200 

8 

8 

0 

1200-2400 

5 

4 

1 

2400-4800 

7 

4 

3 

4800-9600 

6 

4 

6 

AVERAGE  DIFFERENCE' 

2 

13 

“2 - 
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Figure  21.  Drag-Divergence  Curve  for  the  CH-3C 
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SPL  OF  300  -  600  CPS  OCTAVE  MINUS  SPL  OF 
1200-1400  CPS  OCTAVE 
db  (Ref  0. 0002  dyne/tq.  cm,  > 


Figure  23.  Calculated  and  Measured  Levels  -  Hover 


SPL-db(Ref  0. 0002  dyne/sq.  cm.) 
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Figure  24  c.  Calculated  and  Measured  Levels  -  40-Knots  -  Third 

Harmonic 
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Figure  24  d.  Calculated  and  Measured  Levels  -  40-Knots  -  Fourth 

Harmonic 
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SPL-db<Ref  0.0002  <Jyne/*q,  cm.)  ^  SPL-db<R«fO 
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e  25  a.  Calculated  and  Measured  Levels  -  80-Knots  -  First 

Harmonic 
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Figure  25  b.  Calculated  and  Measured  Levels  -  80-Knots  -  Second 

Harmonic 


Figure  25  c.  Calculated  and  Measured  Levels  -  80-Knots  -  Third 

Harmonic 


Figure  25  d.  Calculated  and  Measured  Levels  -  80-Knots  -  Fourth 

Harmonic 
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Figure  26  a.  Calculated  and  Measured  Levels  -  110-Knots  -  First 

Harmonic 


Figure  26  b.  Calculated  and  Measured  Levels  -  110-Knots  -  Second 

Harmonic 


49 


TIP  SPEED  (FT.  PER  SEC.) 

Figure  28.  Calculated  Noise  Levels  at  300  Feet,  20, 000  Founds  Thrust, 

5  Blades 
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Figure  29.  Sound  Pressure  Levels  versus  Thrust  for  3  Tip  Speeds, 

CH-3C 
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OVERALL  VORTEX  ~  (SPL> 


20  24 


28  32 
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THRUST  IN  THOUSANDS  OF  POUNDS 
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Figure  30.  Sound  Pressure  Levels  versus  Thrust  for  3  Tip  Speeds, 

CH-53A 


R*S7«i 


Figure  31*  Trapezoidal  and  Square  Tip  Planforms 


Figure  34.  Comparison  of  Level  at  Strouhal  Frequency  of  Square  and 

Trapezoidal  Tips  @  203  RPM 


Figure  35.  Detection  Level  Criteria  for  Aircraft  -  Spectrum  Level 
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SPL  -  db<R«f  0,0002  dyne/sq,  cm.) 


Figure  36.  An  S-58  Helicopter  Noise  Spectrum  at 

Hover 


ROTATIONAL  NOISE  HARMONIC 

Figure  37.  Effect  of  Blade  Loading  Distribution  on  Noise  Harmonic 

Level 
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Figure  38.  Effect  of  4-per-rev  Blade  Loading  on  Harmonic  Noise  Levels 
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Figure  39.  Effect  of  1-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 
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SPL-db(Ref  0.0002  dync/*q.  cm,) 
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ROTATIONAL  NOISE  HARMONIC 

Figure  40.  Effect  of  2-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 
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ROTATIONAL  NOISE  HARMONIC 


Figure  41.  Effect  of  3-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 
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EOTA-nOHH  iwnumwc 

Figure  42.  Effect  of  4-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 


ROTATIONAL  NO BE  HARMONIC 

Figure  43.  Effect  of  5-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 
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Figure  44.  Effect  of  6-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 


Figure  45.  Effect  of  7-per-rev  Blade  Loading  on  Rotational  Noise 

Harmonics 
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Figure  46.  Influence  of  die  1st  to  the  5th  Blade  Loading  Harmonics 

on  Rotational  Noise  Harmonics 


i 

! 


BLADE  LOADING  HARMONIC 
AMPLITUDE  IS  CONSTANT  AT 


ROTATIONAL  NOISE  HARMONIC 

Figure  47.  Influence  of  the  6th  to  the  10th  Blade  Loading  Harmonics 

on  Rotational  Noise  Harmonics 
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SOUND  PRESSURE  LEVEL-db 


ONLY 

BLADE  LOADING  HARM"'*  I C  ‘*DDED  TO  STEADY 


Figure  49.  Sound  Level  Caused  by  a  Rotor  wu'  .eady  Loading  Plus 
a  Harmonic  Blade  Loading  of  One  Quarter  the  Amplitude  of  the  Steady 
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CONCLUSIONS 


It  as  been  shown  that  single  rotor  overall  vortex  noise  level  for 
square  tipped  blades  with  uniform  inflow  can  be  calculated  with  excellent 
accuracy  by  the  use  of  a  simple  formula  developed  in  this  study.  The 
standard  vortex  noise  spectrum  is  used  for  the  unstalled  condition,  and 
another  spectrum  has  been  developed  for  use  where  stall  is  present. 

This  stall  spectrum  contains  higher  levels  in  the  upper  frequencies.  A 
reduction  in  rotor  tip  speed  or  thrust,  or  an  increase  in  blade  area  (at 
constant  thrust  and  tip  speed)  will  reduce  the  overall  vortex  noise  level. 
Vortex  noise  level  is  highly  dependent  on  tip  shape,  and  substantial 
reductions  may  be  attained  by  proper  design. 

The  solution  derived  for  rotational  noise  may  be  used  to  calculate 
any  number  of  harmonics  in  either  the  near  or  the  far  field  under  uniform 
or  nonuniform  conditions.  Only  single  rotor  systems,  however,  may  be 
considered  using  the  present  solution.  The  only  limitation  to  the  number 
of  harmonics  that  may  be  calculated  appears  to  be  the  availability  of 
accurate  harmonic  blade  loading  data.  Harmonic  blade  loadings  contri¬ 
bute  significantly  to  the  fundamental  rotational  harmonic  and  for  all 
practical  purposes  completely  determine  detectability  levels.  Low- 
frequency  rotational  noise  may  be  reduced  by  lowering  blade  loading  and 
rotor  speed.  High  harmonic  rotational  noise  may  be  reduced  by  reducing 
the  local  stall  and  drag  divergence  tendencies  of  blades  and  by  altering 
aeroelastic  characteristics  to  minimize  harmonic  airloads. 

Blade  slap  is  a  phenomenon  which  affects  both  rotational  and  vortex 
noise  spectra  and  is  generated  by  rapid  changes  of  flow  separation  which 
occur  during  stall  or  drag-divergence  conditions.  Blade  tip  design  is 
extremely  important  in  attenuating  blade  slap.  The  trapezoidal  shaped 
blade  tip  used  in  the  present  study  produced  significant  reductions  in 
blade  slap  noise  level. 
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RECOMMENDATIONS  FOR  FURTHER  STUDY 


1.  Conduct  an  analysis  -  test  correlation  program  to  evalu¬ 
ate  the  accuracy  of  the  rotor  noise  analysis  using  the 
harmonic  airload  results  of  the  NH-3A  (S-61F)  airload 
measurement  program. 

2.  Undertake  further  development  of  existing  aeroelastic- 
aerodynamic  programs  to  extend  their  capabilities  for 
predicting  harmonic  blade  airloads  to  higher  harmonics. 

3.  Perform  a  noise/performance/aeroelastic  trade-off  on  an 
existing  helicopter  to  define  an  optimum  low -noise  con¬ 
figuration.  Modify  a  ship  to  this  configuration  for  verifi¬ 
cation  of  results  and  for  demonstration  purposes. 

4.  Define  the  relationships  between  blade  slap  and  the  degree  of 
stall,  the  extent  of  compressible  drag -divergence,  and  tip 
vortex  characteristics. 
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APPENDIX  I 


SAMPLE  CALCULATION  OF  VORTEX  SOUND  LEVELS 


The  vortex  noise  levels  are  easily  calculated  by  use  of  equation  6 
and  the  vortex  sound  spectra  in  Figure  1. 

Using  the  CH-3C  blade  system  operating  at  203  rpm  and  producing 
thrust  of  19, 000  pounds  for  our  calculations,  the  system  parameters 
listed  in  Table  III  are  used. 


Equation  6  is  used  to  obtain  the  SPL  at  300  feet, 
c  pi  —  I  6.1*10  (^.7)  ■  -  | 

5  r  L  -  io  loy  +2o  100  c  4 

From  Table  III 

Ab=2I7 

Since  for  this  configuration  the  drag  is  small  compared  to  the  lift 
(Cq=0.  077),  it  is  assumed  that  the  lift  is  approximately  equal  to  the 
thrust. 

j  I 

Then  CL=  -  y  (Reference  16) 


Substituting  the  given  values  * 

(n,ooof 

(o.002378)(2I7)(-4«/ 

CL=  0.35 

Substituting  these  values  in  equation  6, 


SPL*  io  Jog  I ♦  20  Jog  0.B7S 


SPL*  81.2  -0.4 
SPL*  80.8  <§>  300 feet 


1 


To  correct  to  distances  of  other  than  300  feet  from  the  rotor,  utilize  the 
expression 

ASPL**olo3^f 

In  this  case  we  wish  to  obtain  the  level  at  225  feet  to  compare  with 
measured  data. 

aSPL*  2olog  fff 


aSPL  =  2.6  db 


Applying  this  correction  to  the  level  at  300  feet,  we  obtain 

SPL- 80.8+2.6  or  S3.4 @  225  feet 

Since  the  operating  range  is  below  stall  on  the  drag-divergence  curve 
of  Figure  21,  the  spectrum  of  Figure  19  should  be  used.  With  a  calcu¬ 
lated  overall  vortex  level  of  83  db  at  225  feet,  the  spectrum  is  then  pre¬ 
dicted  to  be  as  follows: 


Octave  Band 
(CPS) 


SPL 

(db  Ref.  0. 0002  dyne/sq.  cm.) 


150-300  75.0 
300-600  79.0 
600-1200  75.0 
1200-2400  74.5 
2400-4800  70.5 
4800-9600  66.0 


The  measured  values  for  this  condition  were: 


150-300 

75.0 

300-600 

78.0 

600-1200 

74.0 

1200-2400 

74.5 

2400-4800 

70.5 

4800-9600 

66.0 

The  overall  level  for  the  measured  noise  can  be  obtained  by  com¬ 
bining  the  octave  band  levels  logarithmically.  To  obtain  the  overall  SPL 
from  octave  band  levels,  first  determine  the  intensity  ratios  of  the 
individual  octaves.  _ 


SPL. -  10  loq5- 


*Lq  ^0 

Sum  the  intensity  ratios  to  obtain  the  overall  intensity  ratio. 


=  10 


IT.Xi 

VI.  I.' 

Convert  to  overall  SPL. 


SPL  =  io  io9 


The  results  are  summarized  below: 


SPL 

<db) 


Octave  Band 
(CPS) 


Overall 

150-300 

300-600 

600-1200 

1200-2400 

2400-4800 

4800-9600 


Measured 

82.0 

75.0 

78.0 

74.0 

74.5 

70.5 

66.0 


Calculated 

83.0 

75.0 

79.0 

75.0 

74.5 

70.5 

66.0 


APPENDIX  II 


VORTEX  NOISE  TEST  INSTRUMENTATION 


A  block  diagram  of  the  instrumentation  used  in  the  test  program  is 
presented  in  Figure  50.  The  acoustic  signal  was  picked  up  with  a  low 

impedance  moving-coil 


ELECT  iO  -  VOICE 


655C 

MICROPHONE 


REFERENCE  SIGNAL 
GENERAL  RADIO 
TYPE  1307 -A 
TRANSISTOR 
OSCILLATOR 


MEASUREMENT 

GENERAL  RADIO 
TYPE  1551  -C 
SOUND  LEVEL 
METER 


TAPE  RECORDER 


microphone.  The  micro¬ 
phone  was  coupled  to  a 
General  Radio  1551 -C 
Sound  Level  Meter  which 
served  as  a  calibrated  pre¬ 
amplifier  and  decade  attenu¬ 
ator.  The  output  of  the 
sound  level  meter  was 
recorded  on  a  Nagra  III  B 
tape  recorder.  The  over¬ 
all  system  was  calibrated 
electrically  prior  to  each 
series  of  measurements. 


The  recorded  data  was 
played  back  through  an 
octave  band  analyzer.  The 
octave  band  levels  were 
recorded  with  a  graphic 
level  recorder.  Because 
of  the  rapidly  changing 
position  of  the  rotating 
source  and  the  proximity 
Figure  50.  Vortex  Noise  Test  of  the  microphone  with 

Instrumentation  respect  to  the  source,  each 

of  the  octave  bands  showed 
short  time  variations  or  modulation.  In  addition,  there  were  longer  time 
variations  that  were  attributed  to  changing  wind  conditions.  These  time 
variations  were  averaged  out  by  graphically  recording  the  levels  of  the 
individual  octaves  and  noting  the  midpoint  of  the  spread.  The  modulation 
was  found  to  be  higher  at  the  lower  end  of  the  spectrum  than  at  the  high- 
frequency  end  of  the  spectrum. 

The  time-averaged  data  was  then  corrected  for  system  response  and 
the  octave  band  levels  were  tabulated.  The  overall  vortex  noise  level 
was  computed  by  summing  the  corrected  levels  in  the  individual  octaves 
from  the  150-300  cps  band  to  the  4800-9600  cps  band. 


TAPE  RECORDER 


GENERAL- RADIO 
TYPE  1550- A 
OCTAVE  BAND 
ANALYZER 


GENERAL  RADIO 
TYPE  1521 -A 
GRAPHIC  LEVEL 
RECORDER 
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APPENDIX  III 


CALCULATION  OF  PARTIAL  DERIVATIVES 


3S 

Calculate  ~  knowing 

r-2  r~)Z  .  2 


5  -  R  *r  ~  2rR  coscr(cos&co$  T^  +  s  in  &  sin?')  (29) 


From  Figure  3 


>‘°*0=-T&rr  W 


Therefore^  write 


SZ=  x*yi*z1*r*-zrc<>sir,pffipj£srfac<>sr+ysinr) 

1  '  '  (30) 


x*+y*+z* 


which  reduces  to 


-2  2222  ,  ,  .  .  /  \ 

S  -x  +y  +z  +r  -2r(xeosV'  +  ysmir) 

Taking  the  partial  with  respect  to  X 

=x*rcosy' 

but  noting  that  X*  ~RcoS0"cos&  this  reduces  to 
^  cos<r  cosO  ♦  co$y- 

Similarlj^for  the  y  and  H  directions 

^y-*5  ■g"  cos<rsin&~  sir»y' 


(33a) 


(33b) 


(33c) 


77 


It  is  desirable  to  eliminate  the  second  term  of  Equations  (33a)  and 
(33b)  in  order  to  simplify  the  solution.  From  the  ratio  of  magnitudes  of 
the  two  terms  it  appears  that  this  elimination  may  be  valid  if  the  field 
point  distance  ft  is  one  diameter  or  more. 

Equation  (33)  then  becomes 

as  -R  „  a 
37  T  QOi<rcose 

■S’  cost  si’nd 

lz=‘S'  sin<r 
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APPENDIX  IV 


APPROXIMATION  OF  CHORDWISE  PRESSURE  DISTRIBUTION  AS  A 


laajSMNOKBtiMjracKiiiaiasi 


The  Fourier  series  expansion  of  the  initial  blade  passage  time 
history  in  the  rotor  disk  (Figure  5)  is 

AP(r,V,t;=  r  Am  cos  mnnt +  Bw  sinmnjcit 


where 


;  m*l 

t/n+zifan- 

A^- ■***  J Ap(\r}Y",t)  cos  mnnt  4t 

H* 


(34a) 


V/q  4-2101/1 

Bm=^ /a  P(r,V',t)  sin  mn/it  Jt  (34b) 

Wa 

By  shifting  the  time  scale  byH/n  the  pressure  pulse  appears  as  in 
Figure  51.  The  corresponding  mathematical  representation  is 

AP(r,y,t)*  )  am  cos (mn/it- mn'f)*  b^  sin  (mn/it -mn'f) 


*ir/nn 

it/a  P (r/f'#T+  Wa) cos  mnnT  JT 
0 

2i0ui 

bw=  P (r^.T 4*  Vfa)  sin  wn/iT  <1T 


(35a) 


(35b) 


where 

T=t-V/ja 

As  an  approximation  to  the  function  of  Figure  51,  a  step  function 
will  be  used  as  shown  in  Figure  52.  The  amplitude  of  this  step  function  is 

<*+  (36) 
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Figure  51.  Time  History  of  Blade  Passage  Over  rdrdy'After  Time  Shift 


APPENDIX  V 


DEFINITION  OF  BLADE  PITCH  ANGLE 

Blade  angle  f&  is  controlled  by  collective  pitch,  cyclic  pitch,  and 
twist.  With  the  coordinate  system  of  Figure  2,  the  following  is  true: 

longitudinal  cyclic  pitch  angle  =  COS  V" 

lateral  cyclic  pitch  angle  =  $  Sfo* 

where  ft  and  ft  are  the  longitudinal  and  lateral  cyclic  pitch  magni¬ 
tudes,  respectively.  The  effect  of  negative  blade  twist  is  accounted 
for  as  follows: 

pitch  angle  change  due  to  twist  -  -y('--i) 
where 

y  is  the  twist  rate 

r0  is  the  blade  radial  station  at  which  twist  begins. 

The  expression  for  pitch  angle  resulting  from  twist  is  based  on  the 
assumption  that  the  collective  pitch  angle  is  measured  at  the  blade  root. 
Combining  the  effects  of  cyclic  pitch,  twist  effect,  and  the  collec¬ 
tive  pitch  angle,  the  true  blade  pitch  angle  is  obtained. 

p0-T(r-r0)+plcoslf'+ftsm-t- 
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APPENDIX  VI 


DESCRIPTION  OF  PROGRAM  FOR  CALCULATING  ROTATIONAL  NOISE 


A  program  has  been  developed  to  facilitate  the  calculation  of  rotor 
rotational  noise  from  equation  26  of  this  report.  This  is  a  Fortran  IV 
program  intended  for  use  on  the  IBM  7094/44  D.  C.  S.  As  explained 
previously,  the  solution  represents  sound  pressure  at  a  point  for 
rotational  noise  at  the  With  harmonic  of  blade  passage.  The  solution  may 
be  used  for  calculating  the  rms  sound  pressure  at  any  point  in  the  near  or 
far  field  outside  one  diameter  from  the  rotor.  The  solution's  accuracy 
decreases  with  increasing  rotor  system  translational  speed  up  to  Mach 
0. 3  which  is  the  practical  limit  for  accuracy. 

Method  and  Subprograms 

The  method  used  repeatedly  evaluates  the  basic  sound  pressure 
equation  given  in  the  analysis.  This  includes  a  double  integration.  One 
integration  is  around  the  rotor  disk  with  the  sample  points  (azimuth 
angles)  chosen  at  constant  intervals.  Subroutine  SIMCOR,  which  used 
Simpsons  (1/3)  Rule  with  correction  term  was  chosen  to  perform  this 
integration.  The  other  integration  is  along  the  radius  where  sample 
radial  stations  are  unevenly  spaced,  hence,  requiring  the  use  of  sub¬ 
routine  AVQUAD.  This  is  a  method  whereby  several  intervals  are 
approximated  by  "averaged  quadratics"  and  integrated  analytically. 

The  flow  diagram  is  shown  in  Figure  53. 

Running  time  can  be  approximated  at  20  sec.  plus  5  sec.  per  field 
point.  5, 000  units  of  output  is  sufficient  unless  several  causes  are  run 
with  all  output  options  on.  The  deck  setup  is  as  follows: 


$  DCID  with  facility  accounting  information 

$  EXECUTE  IBJOB 

$  IBJOB  MAP 

Binary  Deck  Rotational  Noise 

Binary  Deck  SIMCOR 

Binary  Deck  AVQUAD 

7/8  END  OF  PB 


DATA 
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Figure  53.  Computer  Flow  Diagram 
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7/8  END  OF  PB 
$  IBSYS 
$  PAUSE 
$  PAUSE 


ram  Parameters 


The  parameters  which  must  be  known  to  use  the  program  and  their 
units  are  listed  below: 


n  Number  of  rotor  blades 


Angular  velocity  of  rotor,  revolutions  per  minute 

Blade  chord,  inches 

Velocity  of  sound,  feet  per  second 

Radius  at  start  of  blade  twist,  inches 

Collective  (steady)  pitch  angle,  degrees 

Cosine  component  of  cyclic  pitch  (longitudinal),  degrees 

Sine  component  of  cyclic  pitch  (lateral),  degrees 

Blade  twist  rate,  degrees  per  inch 

Steady  blade  section  loading  for  one  radial  station,  pounds  per 
inch 

Cosine  component  of  the  with  harmonic  of  blade  section  loading 
for  one  radial  station,  pounds  per  inch 

Sine  component  of  the  with  harmonic  of  blade  section  loading  for 
one  radial  station,  pounds  per  inch 

Distance  from  center  of  rotor  to  point  at  which  noise  level  is 
to  be  calculated,  feet 


84 


0  Angle  in  rotor  plane  between  longitudinal  axis  of  helicopter  (0 

degrees  at  tail)  and  point  at  which  level  is  to  be  calculated, 
positive  in  direction  of  rotor  rotation,  degrees 

<r  Angle  in  plane  of  rotor  axis  between  rotor  plane  and  point  at 

which  level  is  to  be  calculated,  positive  in  direction  of  rotor 
thrust,  degrees 

Figures  2  and  3  show  system  parameters. 

Input  Format 

The  input  format,  which  is  to  be  used  for  every  case  is  described 
below.  The  input  format  is  shewn  on  the  coding  form  in  Figure  54.  A 
sample  case  is  shown  in  Figure  55. 

Card  1:  Title  Card 

Col.  1-1,  Col.  2  -  any  desired  title 

Card  2:  Basic  Constants 


Col. 

2  -  end 

code  =  0  last  case;  =  1  case(s)  follow 

Col. 

3  &  4  - 

(right  adjusted)  number  of  field  points 

Col. 

5-12*  - 

n 

Col. 

13-20* 

-  n 

Col. 

21-28* 

-  a 

Col. 

29-36* 

-  c 

Col. 

37-44* 

-R* 

Col. 

45-52* 

-ft, 

Col. 

53-60* 

“ft 

Col. 

61-68* 

A 

Col. 

69-76* 

- r 

*Must  have  decimal  point. 
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ZJTLPX  P*T#r 


otokky  nuNicnrr  shot  t»ui  &a£.  8m*£ 


Card  3:  Control  Constants 

Col.  1  &  2  -  number  of  harmonics  of  rotational  noise 
desired  (l<NS  10  )  (right  adjusted) 

Col.  4  -  number  of  radial  stations  (44  N4  0) 

Col.  5-12  -  angular  increment  of  integration,  degree  (with 
decimal  point) 

(NOTE:  (360/4^)  =  even  integer,  degree 

Card  4:  Output  Controls 

=  00  intermediate  output  not  desired 
=  99  intermediate  output  desired 

Col.  1  &  2  -  controls  printout  of  input  and  radial  integration 

Col.  3  &  4  -  controls  printout  of  azimuthal  integrations 

Col.  5  &  6  -  controls  printout  of  S,/3,u  (from  equation 
26  of  the  report) 

Card  5:  Radial  blade  stations  in  fields  of  9,  inches* 

For  each  radial  station: 

U  in  Col.  1-8*  on  first  card 

Lj  -  L|0  in  fields  of  8*  on  second  card 

M.  -NL  in  fields  of  8*  on  third  card 

I  to 

For  each  field  point: 

One  card  with  three  parameters  in  fields  of  twelve,  R,  B,<r, 
in  that  order  (decimal  points  required). 

*  All  fields  specified  must  contain  points. 
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Notes  Regarding  the  Input  Format 


A  case  (Card  2,  Col.  2)  is  a  complete  set  of  input  data  with  a  number 
of  harmonics,  fields  points,  etc.  A  new  case  would  require  a  new  title, 
a  new  set  of  basic  constants,  etc. 

The  number  of  field  points  here  (Card  2,  Col.  3  &  4)  must  agree 
with  the  number  of  cards  describing  the  individual  field  points  which  are 
the  last  set  of  input  cards. 

The  number  of  rotational  noise  harmonics  desired  (Card  3,  Col.  1 
&  2)  partially  determines  the  program  running  time. 

The  number  of  radial  stations  (Card  3,  Col.  4)  is  arbitrary  and  is 
dependent  on  the  number  of  stations  at  which  blade  loading  is  available. 
Evenly  spaced  stations  over  the  length  of  the  blade  are  desirable  because 
the  program  interpolates  between  points.  When  the  program  is  run  at 
Sikorsky  Aircraft,  the  last  station  is  at  the  blade  tip  and  the  loading  is 
zero  at  this  position. 

The  angular  increment  of  integration  (Card  3,  Col.  5-12)  must  be 
small  to  obtain  a  valid  solution.  If  only  one  harmonic  of  rotational 
noise  is  to  be  calculated,  an  increment  as  large  as  15  degrees  may  be 
used.  When  10  harmonics  of  rotational  noise  are  to  be  calculated^ 
maximum  increment  of  2  degrees  is  tolerable.  Slightly  increased 
accuracy  may  be  gained  by  further  reducing  the  increment,  however, 
computer  running  time  is  increased  accordingly. 

The  intermediate  data  printout  which  is  determined  by  the  output 
control  card  (Card  4)  was  used  primarily  for  debugging  the  program. 

The  normal  printout  includes  the  station  number  and  radius,  azimuth 
number  and  angle,  section  loading  for  each  radial  and  azimuth  position, 
field  point  coordinates,  sound  pressure  and  sound  pressure  level.  The 
intermediate  output  adds  a  great  deal  of  volume  to  the  output  and  is  not 
ordinarily  required. 

The  number  of  radial  blade  stations  (Card  5)  must  correspond  to 
that  specified  in  Card  3,  Col.  4. 

Description  of  Sample  Case 


The  data  for  the  sample  80-knot  flyby  case  was  taken  from  NASA 
TM  X-952,  "A  Tabulation  of  Helicopter  Rotor  Blade  Differential  Pres¬ 
sures,  Stresses,  and  Motions  as  Measured  in  Flight".  The  ship  used 
was  an  H-34  (Sikorsky  S-58)  helicopter.  Rotational  noise  was  calculated 
for  5  harmonics  of  rotational  nou  e  from  the  main  rotor  at  11  field  points. 
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The  blade  root  pitch  ( fa  y/8t ,  {3^  )  was  taken  from  the  pitch  motion  por¬ 
tion  of  Table  13(c)  of  the  referenced  NASA  report.  The  harmonic  loading 
for  7  blade  stations  was  taken  from  Table  13(d)  of  the  same  report  and 
the  tip  was  considered  the  eighth  station  with  zero  loading  assumed.  The 
sample  input  and  output  for  this  case  is  shown  in  Figures  55  and  56, 
respectively. 
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Figure  56  a.  Sample  Output  Sheet 
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Figure  56  b.  Sample  Output  Sheet  (Cont. ) 
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Figure  56  c.  Sample  Output  Sheet  (Cont.) 
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Figure  56  e.  Sample  Output  Sheet  (Cont.) 


APPENDIX  VII 


SAMPLE  CALCULATION-ASSUMED  RELATIONSHIP  BETWEEN  BLADE 


i  M7.<  »j  i  ;  i  mu  n  a 


Referring  to  Figure  7,  which  shows^n  end  view  of  a  rotor  blade, 
and  replacing  cj^with  the  steady  force  F,  note  that 

L -E 

where  ar  is  a  spanwise  element  of  blade.  The  in-plane  force  is  then 

Fj  s  l_Ar  sinfl 

and  the  out-of-plane  force  is 

F^=  Laf  COS/3 

Thrust  and  torque  are  defined  by 

T  =  F^  n  =  nL  Ar  cosfi 
Q  =  F|  rn  =  rnL  Ar  sin£ 

From  the  relationships  we  find 


ta" '  JT 


narcosp 

The  values  selected  for  the  example  are 

0=  7. 79  x  105  pound-inches 
I  =  11, 300  pounds 
r  -  269  inches 
n  =  4  blades 
Ar  -  7  inches 

which  give 


f=  14. 35  degrees 
-  417  pounds  per  inch 
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APPENDIX  VIII 


ROTATIONAL  NOISE  TEST  INSTRUMENTATION 


A  block  diagram  of  the  instrumentation  used  in  the  Phase  II  test  is 
%  presented  in  Figure  53.  Measurements  were  made  with  a  one-half-inch 
condenser  microphone  and  the  signal  fed  into  the  line  input  of  a  Nagra 

I1IB  tape  recorder.  The  400- 
cps  beeper  signal  used  for  a 
time  reference  was  fed  into 
the  recorder  through  its 
microphone  input  and  was 
generated  by  a  General  Radio 
transistorized  oscillator. 
Overall  levels  were  adjusted 
to  avoid  overloading  the 
recorder.  Data  was  recorded 
at  hover,  40-knot,  80-knot, 
and  110 -knot  cruise,  and  70- 
knot  cyclic  pullout.  A  record¬ 
ing  was  made  with  rotor  system 
stopped  and  engine  running  to 
determine  the  extent  of  engine 
noise  interference  in  the  rotor 
noise  frequency  range  of 
interest. 

All  noise  records  were  analyzed  twice.  First,  one-third-octave 
bandwidth  plots  of  noise  level  versus  frequency  were  made  to  determine 
the  frequency  of  the  main  rotor  blade  passage  and  to  evaluate  masking 
noise.  Second,  plots  were  made  of  noise  level  versus  time  at  the  indivi¬ 
dual  rotor  harmonic  frequencies  as  determined  by  the  spectrum  plots. 

Frequency  analysis  of  the  records  with  the  engine  running  and 
stopped  rotor  system  showed  that  engine  spikes  appeared  from  79  cps 
upward.  The  tail  rotor  fundamental  blade  passage  frequency  of  the  S-58 
helicopter  appears  at  85  cps  for  the  rotor  speed  used  in  the  test.  The 
fifth  harmonic  of  the  main  rotor  is  at  75  cps  and  is  masked  by  engine 
noise.  For  this  reason,  wily  the  first  four  harmonics  of  the  rotational 
noise  were  analyzed.  This  data  is  presented  in  the  section  entitled 
"Correlation  of  Measured  and  Calculated  Rotational  Noise". 


GENERAL  RACK) 
t m  lisi  fiL 

COMDEME*  MlCRGfHQMB 
SYSTEM 


ANALYSES 


Figure  57.  Rotational  Noise  Test 
Instrumentation 
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APPENDIX  IX 


ROTOR  SYSTEM  PARAMETERS  AND  LOADS  FOR  PHASE  III  STUDY 


The  S-58  rotor  system  in  the  hover  condition  is  used  as  the  basis 
for  the  Phase  III  study  of  the  effects  of  harmonic  blade  loading.  The 
pertinent  rotor  system  parameters  are: 


Number  of  blades  4 

Blade  radius,  feet  28 

Blade  twist,  degrees  -8 

Blade  chord,  feet  1.367 

Selected  operating  parameters  are: 

Rotor  speed,  rpm  210 

Blade  pitch  angle  at  root,  degrees  19. 6 

Thrust,  pounds  11,300 

Loaded  radial  increment,  inches  7.0 


Eight-tenths  of  the  tip  radius  (r  =  269  inches)  is  selected  as  the  effective 
blade  radius.  The  pitch  angle  £ ,  corrected  from  the  root  angle  for 
blade  twist  is  14.35  degrees. 


From  Figure  58, 

L  =  blade  section  loading,  pounds 
per  inch 

r  *  loaded  blade  incremental  radius, 
inches 

P  =  blade  pitch  angle,  degrees 
-  number  of  blades 

Total  thrust ”T  =  «L  Ar  C0S0 
Figure  58.  Blade  Aerodynamic  . 

Loading  or,  solving  for  L, 

l  |  _  11.303 

L*  ntrcosp  *  4*7«C0SH.JJ*  =  417  pounds  per  in. 

This  417  pounds -per -inch  blade  loading  is  used  as  the  steady  component 
of  all  Phase  III  calculations.  Harmonic  loading  levels  are  fractions  of 
this  figure. 
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